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Abstract 

We present a new method for solving symbolically zero-dimensional polynomial 
equation systems in the affine and toric case. The main feature of our method is the use 
of problem adapted data structures : arithmetic networks and straight-line programs. 
For sequential time complexity measured by network size we obtain the following 
result : it is possible to solve any affine or toric zero-dimensional equation system in 
non-uniform sequential time which is polynomial in the length of the input description 
and the "geometric degree" of the equation system. Here, the input is thought to be 
given by a straight-line program (or alternatively in sparse representation), and the 
length of the input is measured by number of variables, degree of equations and size of 
the program (or sparsity of the equations). The geometric degree of the input system 
has to be adequately defined. It is always bounded by the algebraic-combinatoric 
"Bezout number" of the system which is given by the Hilbert function of a suitable 
homogeneous ideal. However, in many important cases, the value of the geometric 
degree of the system is much smaller than its Bezout number since this geometric 
degree does not take into account multiplicities or degrees of extraneous components 
(which may appear at infinity in the affine case or may be contained in some coordinate 
' hyperplane in the toric case). 

Our method contains a new application of a classic tool to symbolic computa- 
tion : we use Newton iteration in order to simplify straight-line programs occurring 
in elimination procedures. Our new technique allows for practical implementations a 
meaningful characterization of the intrinsic algebraic complexity of typic elimination 
problems and reduces the still unanswered question of their intrinsic bit complexity 
to algorithmic arithmetics. However our algorithms are not rational anymore as are 
the usual ones in elimination theory. They require some restricted computing with 
algebraic numbers. This is due to its numeric ingredients (Newton iteration). Never- 
theless, at least in the case of polynomial equation systems depending on parameters, 
the practical advantage of our method with respect to more traditional ones in sym- 
bolic and numeric computation is clearly visible. Our approach produces immediately 
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a series of division theorems (effective Nullstellensatze) with new and more differenti- 
ated degree and complexity bounds (we shall state two of them). 

It should be well understood that our method does not improve the well known 
worst-case complexity bounds for zero-dimensional equation solving in symbolic and 
numeric computing. 



Part of the results of this paper were announced in 
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1 Introduction and Statement of Results 

Let k be an infinite, perfect field which we think "effective" with respect to arithmetic 
operations as addition/subtraction, multiplication/division and extraction of p— th roots 
in case that k has positive characteristic p . Let A; be a fixed algebraic closure of k . 

An important problem in elimination theory is the computation of the isolated points 
of an affine algebraic variety. This means we consider the algebraic variety 

V:=V(ft,..., f s ) := {x G A n (k) : f x {x) = 0, . . . , f s (x) = 0}, 

where fx, . . . ,f s G k[X±, . . . , X n ] with s > n are polynomials of degree at most d in the 
variables X\ , . . . , X n and we look at the following problem : given a nonzero linear form 
H G k[X\, . . . , X n ] and a new variable T , compute a nonzero polynomial p £ k[T] such 
that p(H) vanishes on all isolated points of V . This "weak" form of symbolic solving 
of the equation system given by the polynomials fx , . . . , f s is in all thinkable aspects 
(effectivity, practical efficiency and theoretical complexity) equivalent to the following 
more explicit form : find univariate polynomials q, v i, . . . , v n G k[T] such that the set of 
isolated points of V can be written as {(vx(t), v n (t)) ; t G A 1 (A;) , q(t) = 0} (see 
[24|, @ and Lemma [7| below). 

As many authors do (see for instance, §, |], §, §, @, @, 0, 0], @, @, 



p9fl , p0[ , M], p3[ |, plj ], [53]), we replace the original input system by n generic 
k -linear combinations of the equations fx , • • • , f s ■ This preparation of the input con- 
serves all irreducible components (and, in particular, the isolated points) of V and adds 
possibly some new, extraneous ones (these components have to be eliminated afterwards). 
Thus, let us suppose from now on that s = n holds. This implies that the input sys- 
tem fx,...,f n forms locally a complete intersection with respect to the isolated points 
of V . In view of the surprisingly low (linear) regularity bound for homogeneous global 
complete intersection ideals (see || , j|4| , |l6| ) one is tempted to replace the original equa- 
tions fx, ■ ■ ■ , f n G k[Xx, . . . , X n ] by new homogeneous ones, G\, . . . , G n G k[Xo, . . . , X n ] , 
and the affine variety V = V(fx, ■ ■ ■ , f n ) by the corresponding projective one, namely 
W := V(Gx,...,G n ) = {x G TP n {k) ; Gx{x) = 0, . . . , G n (x) = 0} . A simple minded 
way to proceed consists in taking as G%, . . . , G n just the homogenizations of the original 
system fx , ■ ■ ■ , f n ■ However this idea does not produce the desired effect even in case that 
the polynomials fx, ■ ■ ■ , fn form a regular sequence in k[Xx, ■ ■ ■ , X n ] . The reason is that 
the homogenizations of fx, ■ ■ ■ , fn need not to generate a complete intersection ideal. 

In order to remedy this defect, many authors use homotopic deformations and path- 
following methods for the construction of a regular sequence of auxiliar homogeneous 
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polynomials G\, . . . , G n , from whose zeroes the solutions of the original system fi,...,f n 



are then extracted (see [|, 0, Q, @, U, ||] ; [24|, [27|, EI, Q, @, [47|, 



||5q| ). The polynomials Gx, . . . ,G n are homogeneous in the variables Xq,Xx, . . . , X n . 
Moreover they depend on a deformation parameter e and form a regular sequence in 
k(e)[Xo, Xi, . . . ,X n ] . They define a zero dimensional projective subvariety of TP n (k(e)) 
without points at infinity (here k(e) denotes any algebraic closure of the function field 
k(s) ). Specializing the polynomials Gi, . . . , G n in e = we obtain a projective subvariety 
of TP n (k) which contains the isolated points of V as irreducible components. Since the 
underlying deformation is flat in the isolated points of V , well known techniques based 
on the implicit or explicit use of Macaulay's w-resultant allow first to find the solutions of 
V(G±, . . . , G n ) and finally the isolated points of V . 

This deformation method introduces a somehow artificial dependency of the complex- 
ity of the algorithms on the regularity of the Hilbert function of the homogeneous ideal 
(Gi, . . . , G n ) of k(e)[Xo, X\, . . . , X n ] and on the degree of (G±, . . . , G n ) which is defined 
by means of the Hilbert polynomial and is called its Bezout number. The regularity of 
(G\, . . . , G n ) is bounded by nd — n. This implies that deformation based algorithms have 

to triangulate rectangular matrices of size ( ] x n ( ] in order to extract from this 

data some "geometrically meaningful" square matrix of size the Bezout number. The 
characteristic polynomial of this square matrix is the basic eliminating form that all al- 
gorithms work with. However, the Bezout number of the ideal (Gi, . . . , G n ) (i.e. the 
degree of the eliminating form) includes multiplicities and degrees of extraneous compo- 
nents which the previous deformation process adds to the original variety V . The exact 

n 

value of the Bezout number is Yl degGi which is of order d n . 

i=l 

The method we present in this paper will be independent of "algebraic" quantities such 
as regularity of the Hilbert function or Bezout number of an appropriate homogeneous 
ideal. This allows in geometrically well suited cases (typically when the geometric degree 
of the variety V is low) to reduce the size of the matrices in the algorithms (and hence the 
complexity of the procedures). However, in worst case (e.g. when the equations fi, ■ ■ ■ , f n 
are generic) our complexity bounds are roughly the same as those obtained by deformation 
based algorithms. This is not surprising in view of the known lower bounds for elimination 
problems ( Q, §7|]). 

An important aspect of elimination procedures consists in the encoding of the poly- 
nomials appearing as inputs, outputs or intermediate results of the algorithm. Encoding 
polynomials by their coefficients (dense representation) faces us with an input size of order 

(d _|_ n \ 
< c • d n (if n, d > 3 we can take c = 1 , otherwise we can take as c the base 
n J 

of the natural logarithm). If the equations fi, . . . , f n are generic, the geometric degree 
of the variety V and the Bezout number of the system (given by the homogenizations of 
fx, . . . , f n ) coincide and are of order d n . Thus, a complexity bound of type d°^ is both 
polynomial in the dense input size and the degree of the variety V , what is the best we 
can hope in worst case (see jn]] and |57]]). 

However, when the geometric degree of V is low, one wishes to use more economic 
encodings of inputs, outputs and intermediate results. Such an encoding is given, for 
instance, by the data structure straight-line program. 
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For the moment let us fix the following notations and assumptions : 
There is given a family of n input polynomials fi, ■ ■ ■ , f n £ k[Xi, • • • > -^n] which are 
thought to be encoded by a straight-line program (arithmetic circuit) (3 without essential 
divisions in k\X\, . . . ,X n ] (this means we allow [3 only to contain divisions by non-zero 
elements of k). We denote the nonscalar size ("nonscalar length" in more traditional 



terminology) and the nonscalar depth of (3 by L and I respectively (see [23], [30], [40], 
[59] and Section |2| for the notions of straight-line program, arithmetic network and the 
complexity measures we shall use subsequently). 

Let us first consider the affine case. Here, we suppose that fi, ■ ■ ■ , f n form a regular 
sequence in k[X\, . . . ,X n ] . Fix 1 < i < n. The affine variety defined by the ideal 
(/i, namely V(/ij •••>/*) > is a Zariski closed subset of A n := A n (k) of pure 

dimension n — i (i.e. all irreducible components of V(fi, ■ ■ ■ , fi) have the same dimension 
n — i; we say also that V(/i, • • • , fi) is an equidimensional variety of dimension n — i). 
The (geometric) affine degree of V(/i, ...,/») is defined as usual as the cardinality of the 
finite set of points we obtain cutting V(f%, . . . , fi) by n — i generic affine hyperplanes of 
A n (more generally, we define the affine degree of a closed Zariski subset of A™ as the sum 
of the degrees of its irreducible components. See e.g. |29j or for this notion of degree 
and its motivation). We denote the affine degree of V(/i, • ■ ■ , fi) by degV(f\, . . . , fi) . 

Let us now define 5 := max{degV(fi, . . . , fi) ; 1 < i < n} as the (geometric) affine 
degree of the equation system f\, ...,/„ . We write V := V(fi, . . . , f n ) and observe that 
V contains only finitely many, namely degV < S points. 

Let us now consider the toric case. Fix 1 < i < n . We define the toric irreducible 
components of the affine variety V(/i, . . . , fi) as those which are not contained in one 
of the hyperplanes V(Xj) , 1 < j < n , of A n . We suppose that fx, ■ ■ ■ , fi form a toric 
complete intersection. This means that toric irreducible components of V(f\, . . . , fi) exist 
having all dimension n — i . The (geometric) toric degree of V(/i, • • • ,fi) is denoted by 
deg*V(fi, . . . , fi) and is defined as the sum of the affine degrees of all toric components 
of V(/i, . . . , fi) . We call 5* := max{deg*V(fi, . . . , fi) ; 1 < i < n} the (geometric) toric 
degree of the equation system f\ , . . . , f n . 

Let be given a non-zero linear form H of k[Xi, . . . ,X n ] represented by its coefficients 
and let T be a new variable. We consider as the fixed input of all our principal algo- 
rithms the division free straight-line program (3 which computes fi,...,f n and, where it 
makes sense, the coefficient n -tuple of H . We recall that d is an upper bound for the 
degrees of the polynomials of the input system fi, . . . , f n , that 5 and 5* represent the ge- 
ometric affine degree and the geometric toric degree of the system fi, ■ ■ ■ , f n respectively. 
Furthermore we recall that L and £ are the nonscalar size and depth of (3 respectively. 

With these notations and assumptions, we may state our main results as follows : 



Theorem 1 (the affine case) Assume that fi, ■ ■ ■ ,f n form a regular sequence in 
k[X\, . . . , X n ] . Suppose furthermore that for 1 < i < n the ideal (/i, . . . , fi) generated 
by fi,...,fi in k[X\, . . . , X n ] is radical. Then there exists an arithmetic network with 
parameters in k which, from the input given by (3 and H , computes the coefficients of a 
monic polynomial p G k[T] such that p(H) vanishes on V = V(/i, • • • , fn) ■ The nonscalar 
size and depth of the network are (ndSL)°^ and 0(n(log2(nd) + l)log2$) respectively. 
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Theorem 2 (the toric case) Assume that for 1 < i < n the polynomials /i, •••,/« 
form a toric complete intersection and generate an ideal in k[X\, . . . ,X n ] whose localiza- 
tion (fx, . . . , fi) n is radical. Then there exists an arithmetic network with parameters 

3=1 

in k which, from the input given by (3 and H , computes the coefficients of a monic poly- 

n 

nomial p* G k[T] such that p*(H) vanishes on the toric variety V* := V \ V( JJ Xj) . 

j=l 

The nonscalar size and depth of the network are (nd8* L)°W and O(n(log2(nd)+£)log20~*) 
respectively. 

Let us remark that the algorithms underlying Theorems [l] and [2] make substantial use 
of linear algebra subroutines dealing with matrices of size at most 2dS (affine case) and 
2d5* (toric case). 

We may sharpen the assertions of Theorem [l] and ^ to the following statement : 

Theorem 3 Let fi,...,f n be polynomials of k[X\, . . . , X n ] which satisfy the assump- 
tions of Theorems [7] or [^. Then there exists an arithmetic network with parameters in k 
which from the input given by the circuit (3 computes in the affine case the coefficients 
of polynomials q, v\, . . . , v n S k[T] and in the toric case the coefficients of polynomials 
q*,v*, . . . ,1;* S k[T] such that the following conditions hold : 

i) q and q* are monic and have degrees degq = degV and degq* = degV* = 
deg* V respectively. Moreover vi,...,v n and v[* , . . . , satisfy the degree bounds 
max{degvi ; 1 < i < n} < degq and max{degv* ; 1 < i < n} < degq* . 

ii) the zerodimensional varieties V and V* can be parameterized as 

V= {{ Vl (t),...,v n {t)) ; t G A 1 , q(t) = 0} and 
V*= {(t>i(t), • • • ,w*(i)) ; t £ A 1 , q*(t) = 0} 

The nonscalar size and depth of the network are (nd5L)°^ and 0(n(log2(nd) +£)log2S) 
in the affine case and (nd5* L)°^ and 0(n(log2(nd) + £)log20~*) in the toric case. 

Of course we can reformulate Theorems [l], ^ and || for sparse inputs. In this case 
the parameters L and t measuring nonscalar size and depth of the input circuit (5 
can be replaced by the sparsity N of the polynomials /i,...,/ n and by log2d respec- 
tively (the sparsity is the number of monomials with non-zero coefficients appearing in 
fi, . . . , f n ). The nonscalar size and depth of the resulting network are then (ndSN) ^ 1 ' 
and O(nlog2(nd)log20~) in the affine case and (ndS*N)°^ and 0(nlog 2 (nd)log 2 6*) in the 
toric case 

The most important context of applications of Theorems [j], § and || is the following 
situation of "parametric/numeric equation solving" : 

let Q be a prime field with algebraic closure Ct and let 9i,...,9 m be indetermi- 
nates over (called "parameters"). Let k := 0(#x, • • • ,9 m ) (there will be no harm to 
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our arguments in characteristic p since Q is perfect). Suppose now that (3 is a divi- 
sion free straight-line program in Q[9i, . . . , 9 m , Xi, . . . , X n ] of nonscalar length and depth 
L and £ respectively. Thus our polynomials fi, ■ ■ ■ , f n belong to the polynomial ring 
. . . , 9 m , X\, . . . , X n ] . Assume that H is a nonzero linear form of tl[Xi, . . . , X n ] . 
The varieties V and V* are interpreted as subvarieties of the affine space A n (k) = 
A™(f2(#i, . . . ,9 m )) ■ The remaining assumptions and notations are the same as in Theo- 
rems %%and^for k := Q(0i, . . . ,9 m ) . 

In this situation the statement of Theorem [I], [2] and ^ can be sharpened as follows : 

Theorem 4 There exists an arithmetic network with parameters from the field , which 
has nonscalar length (nd5L)°^ and depth 0(n(log2(nd) + £)log20~) in the affine and 
which has nonscalar length (nd5*L)°^ and depth 0(n(log2(nd) + £)log2d~*) in the toric 
case. This arithmetic network produces as output a division free straight-line program in 
Cl[9i,. . . ,9 m ] , say 7, in the affine and 7* in the toric case. This straight-line program 
has the following properties : 

• according to the case (affine or toric) the circuits 7 and 7* represent the coeffi- 
cients with respect to the variable T of polynomials p and p* or q,wi, . . . ,w n and 
q* ,w*, . . . , w n which belong to . . . , 9 m , T] . Moreover in the situation where it 
makes sense 7 and 7* compute also the (nonzero) discriminants p G . . . , 9 rn \ 
and p* E Q[9\, . . . , 9 m ] of q and q* . 

• the straightAine programs 7 and 7* have asymptotically the same nonscalar size 
and depth as the arithmetic network generating them. 

• the (monic in T ) polynomials p,p*,q,q* £ Q[9±, . . . ,9 m ,T] and the rational func- 
tions Vl :=f € fi(0i, . . . ,9 m )[T] and vf := 2* . . . , 2& € fi(0 x , . . . , 9 m )[T] 
have the properties stated for them in Theorems [J, and [| respectively. For instance 
in the situation of T heorem |j| the polynomial p{6\, . . . , 9 m , H) vanishes on V = 
V(fi, . . . , f n ) = {x £ f2(0i, . . . , 9m) ; f\(x) = 0, . . . , f n {x) = 0} and in the situa- 
tion of Theorem^(i) the algebraic variety V has the form V = {v\(t), . . . ,v n (t) ; t £ 

,6 m ) , q(9 l ,...,9 m ,t)=0}. 

We shall omit the proof of Theorem |I| because it would be almost textually the same 
as that of Proposition |2| below. 

A natural question to ask is the following one : what "real life" complexities are hidden 
behind the notions of arithmetic network and straight -line program with parameters in 
k? What is the link of these data structures to "ordinary" rational arithmetic networks 
and straight-line programs with parameters in the coefficient field k? For answering 
this question let us restrict ourselves to the affine case (the toric case can be discussed 
analogously) . 

A simple minded translation of the arithmetic network with parameters in k underly- 
ing Theorem [l] to the rational context over k would produce an ordinary arithmetic net- 
work with parameters in k of nonscalar size and depth (nd5 n L)°^ and 0(n(log 2 (d5)+£)) 
respectively. Whereas the nonscalar depth is fine, the nonscalar size of the network grows 
up to order of the sequential time complexity of usual Groebner basis computations for 
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the elimination task we are considering, i.e. to {nd n )°« in worst case (here we assume 
that 5 = d n and L = ^ j = 0{d n ) holds; see |] and @). 

The problem arising in this way can be localized in the iterative character of the algo- 
rithm (n iterations) which in the ordinary rational arithmetic network version produces 
stupidly growing straight-line programs for the representation of intermediate results. 
This is due to the repeated use of interpolation subroutines which nevertheless can be 
summarized by suitable FOR commands (see [25]). However, applying "data compres- 
sion" by means of the already mentioned homotopic deformation method one can lower 
this ordinary network size to (nd n )°W which is an already known complexity bound for 
the elimination task we are considering (@, Q, ||], @, @. See also @, @, g|, 
[44]). This means the advantage of the rational version of our method is reduced to ex- 
egetic "practical complexity" as far as running time is concerned (nevertheless we may 
have saved something with respect to storage space). 

However a more radical compression of the enormously expanding straight -line pro- 
grams during the execution of the algorithms is possible thanks to a symbolic adaptation 
of Newton's method (or Hensel iteration as someones like to call it in case of positive char- 
acteristic) in combination with "Vermeidung von Divisionen" (|30|). This compression is 
done by the algorithm underlying the fundamental technical Lemmas |l] and § below. 

The key point of our method consists in a consequent use of problem adapted data 
structures : arithmetic networks and straight-line programs with parameters in k (or 
k). A similar ground idea may be found in the impressive theoretical and practical work 
of D. Duval and her school on dynamic evaluation (fllll). Closest to our complexity 
results comes the numeric method of M. Shub and S. Smale for finding "approximate" 
solutions of zero-dimensional homogeneous "average" equations systems (]56]]). Their 
method too is independent on the Bezout number and the regularity of the input ideal. 
The sequential time complexity of their algorithm (measured as ours by counting just 

(the size of the densely given input 

n J 

system) and the number of approximate solutions wanted. However caution is necessary 
in the interpretation of their result as "equation solving" in the sense of computer algebra. 
First they need sufficient genericity of the system (just the contrary of the point of view 
of computer algebra which focuses on "special systems"). Secondly Newton iteration (for 
which they seek approximate solutions) is inefficient in terms of bit-complexity if the 
rational numbers appearing in it are given in binary coding. This fact is due to Liouville's 
estimate (see J50[|). However this second drawback can be avoided coding rational numbers 
by straight-line programs. Identity of rational numbers given by such encodings can be 
detected in random polynomial time by a BPP-test (PH, p4|]) or a fixed sample test (as 



in 1 32] or |L3[]) which is due to [pq| . Summarizing this we can say that in case $7 := Q 
and k:= Q(#i, . . . , 9 m ) with 6\, . . . , 6 m algebraically independent, our algorithms have a 
reasonable translation to the (probabilistic) bit-model and that this translation conserves 
the overall complexity character of our algorithms (this is work in progress and will be the 
subject of a further publication; see Section [?]). Of course it would also be advantageous 
to dispose of algorithms of the same complexity type which are realizable exclusively by 
rational arithmetic networks and straight -line programs (with parameters in k). This is 
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a question which we shall also consider in further work. 

Our paper is organized as follows : 
in Section [2] we introduce our computational model of arithmetic networks and straight- 
line programs with parameters in k (or k ). In Section ||| we prove Lemma |l| and § which 
represent a key tool for our algorithms. Thanks to these lemmas we are able to compress 
efficiently certain straight-line programs which appear as intermediate computations in 
our procedures. The algorithm underlying Lemma [l] and ^ requires the use of algebraic 
numbers (observe that this situation produces the side effect that our subsequent arith- 
metic networks and straight-line programs depend on algebraic parameters). In Section 
we review some known technical lemmas which are necessary for the proofs in Section 
The essential steps of the proof of Theorem [l], ^ and [3| are contained in Section [|. 
This proof is based on a recursive construction of a Noether normalization and a primitive 
element for the varieties V(f\, . . . , fi) , 1 < i < n, followed by a procedure of cleaning 
extraneous irreducible components. This recursive construction is explained in Section || 
and a recursion law of the complexity of the underlying algorithm is given. At this point 
the application of Lemma [l] and [2] of Section ||| becomes crucial. A simple minded iteration 
of the recursive construction we introduce in Section || would lead to a size explosion of 
straight -line programs. Lemma [l] and |2] allow us to compress after each iteration step the 
straight-line programs occuring in our algorithm. The core of the proof of Theorem |l|, 
H and H is contained in Proposition |2] and ||] which generalizes the statement of Theorem 
^ to the case of arbitrary reduced complete intersection ideals. This proposition is used 
in Section || in order to formulate two division theorems (effective Nullstellensatze) with 
new complexity and degree bounds which are polynomial in our parameters n, d, 5 and L 
(Theorem |5| and |6|). These new division theorems do not improve the degree (and height) 
bounds of the classical ones ([|, §, g, §, 0, @, @, fyj, H) but they have 



three advantages : the first one is that they are better suited for computational issues, the 
second one is that they are more general and the third one is that they explain better the 
classical results. 



2 The Computational Model 



The idea of using straight -line programs as succinct encodings of special polynomials (as 
they appear in elimination problems) goes back to the late seventies and was discovered 
by different people independently (one of them being our coauthor J. Morgenstern who 
prefered oral "hadise" to hardcover publications for the dissemination of his thought). 
This idea appears implicitly or explicitly in the following representative (although not ex- 
haustive) list of papers, mostly dedicated to probabilistic testing of polynomial identities : 
p2[ , |34|| , J55[ , [15], [64]. First applications of this point of view to computer algebra dealt 
only with elimination of just one variable (absolute primality testing, greatest common 
divisor computation and factorization of multivariate polynomials; see |3J 



|3|, §§) and 

were later extended to multivariate elimination problems by means of "duality techniques" 
of different type and "Vermeidung von Divisionen" (@, §E§, §§, @, gp; see also |l|). 



Straight -line programs and arithmetic circuits have their origins in the design of semi- 
numerical algorithms and represent a link between numerical analysis and computer sci- 



ence (see [23], [[}(]], [59] and the references given there). In this section, we slightly modify 
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the notion of "ordinary" arithmetic networks adding a special type of computation nodes, 
called "algebraic gates" . These gates display (generally algebraic) elements of k which are 
given as zeroes of a precomputed univariate polynomial over k . We are now going to ex- 
plain what we mean by "arithmetic network" and "straight-line program with parameters 
in k (or k)" . Let for the moment K be any of the fields k or k . 

A piecewise rational function (over K ) is a mapping (p : A n — ► A 1 such that there 
exists a partition of A n in K -definable constructible subsets {Cj : 1 < j < M} , such 
that for any 1 < j < M there exists a rational function <j)j £ K(X\, . . . , X n ) defined 
everywhere on Cj , which verifies : 



f I r = Pj I n ■ 

An " ordinary" arithmetic network over A' is a device that evaluates piecewise rational 
functions. However, in our applications, most of the functions are piecewise polynomial, 
i.e. the <pj are polynomials belonging to K[Xi, . . . , X n ] . 

An arithmetic network T over A is a pair T = (G,Q), where Q is directed acyclic 
graph, with vertices (called nodes or gates of Y) of indegree 0, 1, 2 or 3 and where Q is 
an assignment (labeling) of instructions and (piecewise) rational functions ( "intermediate 
results") to nodes (the labeling will be specified subsequently). The graph Q contains n 
nodes of indegree 0, labeled by the variables Xi,..., X n , which are called the input nodes 
of r. 

We define the depth of a node v of the graph Q as the length of the longest path 
joining v with some input gate. Let us denote any node of Q by a pair of integers (i, j) , 
where i represents the depth of the node and j is an ordinal number assigned to the node 
which is given by some ordering of the set of nodes of depth i (see [4C], [48] for this type 
of encoding of arithmetic networks). 

Let be a node of Q . If has indegree and does not represent an input, 

the labeling Q assigns to this node a constant from K (which we call subsequently a 
parameter of V). If has indegree 2 the labeling Q assigns to it an arithmetic 

operation of k. If has indegree 1 the labeling Q assigns to it a sign (i.e. equal to 

zero) test and a boolean variable. If (i,j) has indegree 3 the labeling Q assigns to it a 
selector which makes a choice between two prefixed nodes of Q according to the value of 
a boolean variable associated with a third node. A node of Q labeled in this manner is 
called a gate of T . We have already introduced the input gates. A gate which is labeled by 
a constant or an arithmetic operation of K is called a computation gate (in case of positive 
characteristic p we consider also extraction of p-th. roots of elements of K as arithmetic 
operations). The other gates are called sign or selector gates, depending on their labeling. 
We suppose that the labeling Q is "meaningful", i.e. computation gates have ingoing 
edges (if any) coming from computation or input gates, sign gates have an ingoing edge 
coming from a computation gate , and selector gates have two ingoing edges coming from 
computation or input gates and another ingoing edge coming from a sign gate. Under this 
hypothesis the labeling Q associates in a obvious way to each computation gate (i,j) a 
rational function which we denote by Qij . Let (i, j) be a sign gate with ingoing edge 
coming from a node (r, s) and let x be a point of A™ where Q T)S is defined. Then the 
boolean variable Bij associated to (i,j) takes the value 1 or according to the truth or 
falseness of the statement "Qij(x) =0". If (i, j) is a selector gate with two ingoing edges 
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coming from computation gates (k,l) and (k',l') and if x is a point of A n where Qki 
and Qk' v are defined and if the third ingoing edge of (i, j) comes from the sign gate (r, s) 
then we associate to the field element Qk,i{x) or Qk> v(%) according to the value 

which takes in x the boolean variable B TyS associated to the sign gate. We consider the 
gates of r with outdegree as outputs and we suppose that for each x € A n there exists 
a "stream" of "consistent" paths (with all rational functions defined in x ) to a prefixed 
number of output gates. 

A computation gate of V is called nonscalar (with respect to ) if the instruction 
associated to it corresponds to a multiplication of nonconstant rational functions or to a 
division by a nonconstant rational function. We assign to nonscalar computation gates 
unit costs whereas all other gates are taken cost-free (in particular, computation gates of T 
corresponding to K -linear operations are free). In this way, we associate two complexity 
measures to the arithmetic network T : 



nonscalar sequential time or the nonscalar size of V , defined as the total number of 
nonscalar computation gates of T , 

nonscalar parallel time or the nonscalar depth of V defined as the longest oriented 
path of r joining an input gate with an output gate when only nonscalar computa- 
tion gates are taken into account for the "length". 

These two nonscalar complexity measures are fairly realistic with respect to operation 
counting because V can always be rearranged such that the total number of gates and 
the total depth of T are bounded roughly by the square of its nonscalar size and depth 
respectively. The intuitive meaning of the nonscalar size of V is (sequential) running time 
whereas the nonscalar depth is linked to rather mathematical quantities as degree and 
height (in case K := Q) of the rational functions appearing as intermediate results of T . 



The total depth of T has a natural interpretation as minimal storage space (see [4C] for 
more details on this subject). 

A special case of arithmetic networks are those which have as single output a boolean 
combination of sign gates. We call them decision networks. Another important case of 
special networks is represented by those which contain only computation gates (no sign 
or selector gates). They are called straight-line programs or arithmetic circuits and they 
compute (or represent) rational functions belonging to K(X\, . . . , X n ) . Often we call the 
(nonscalar) size of a straight-line program its (nonscalar) length. 

An arithmetic network or a straight -line program without any nonscalar division (only 
divisions by non-zero elements of K are allowed in this case) is called division free. 
Division free straight -line programs compute polynomials of K[X\, . . . ,X n ] . Finally we 
say that a family of arithmetic networks with nonscalar sequential time cost function L is 
well parallelizable if its nonscalar parallel time cost function I behaves as I = 0(log<it>) . 

Since we are not considering other complexity measures than the nonscalar ones, we 
simplify our terminology from now on saying just size/ sequential time or depth /parallel 
time with reference to the nonscalar cost measure. 

Suppose now K := k . We are going to extend the model of "ordinary" arithmetic 
networks (with parameters in k ) introducing a new type of computation nodes, called 
algebraic gates. These algebraic gates involve for given inputs from k elements of k which 
generally are algebraic over k ( in other words, these gates "compute" algebraic functions 
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in the inputs). In order to explain the nature of these new gates let us assume that the 
arithmetic network T in question disposes over a second type of input nodes, labeled by 
indeterminates, say A±, . . . , A m which are called parameter variables. To understand this, 
remember that the concrete input for our geometric problems (e.g. in Theorem |], [2] and 
^) is always a division free straight -line program (3 in k[Xi, . . . , X n ] which represents the 
equations of the input system. The parameter variables Ai,... ,A m of T are specialized 
for a concrete input circuit (3 into the values of the program parameters of (3 , i.e. into 
values belonging to k . Therefore in all our applications the number m will be of order 
L 2 . Some computation nodes of the arithmetic circuit T will therefore depend exclusively 
on the parameter variables A%, . . . , A m . The labeling Q assigns to these nodes a rational 
function of k(A\, • • • , A m ) . We call such a node of T parameter gate. Let T be a new 
indeterminate. An algebraic gate of T is now a node (i,j) which has indegree and 
outdegree N , where N is an arbitrary (but fixed) natural number. The ingoing edges of 
(i, j) are supposed to come from TV" parameter gates of T , say (sq, ro), . . . , (sjy-i, "nv-i) , 
to which the labeling Q assigns polynomials /io, . . . , h^-i £ k[A\, . . . , A m ] . Consider now 
an arbitrary specialization a = (oci, . . . , a m ) £ k m of the parameter variables A%, . . . , A m 
in k (such a specialization represents in our applications always a concrete input circuit 
f3). The polynomials ho, ... , /ijv_i take in the argument a values r/o := ho(a), . . . , r]N-i ■= 
hN-\{a) which belong to k. The algebraic gate (i,j) assigns then to its N outgoing edges 
the totality of the zeroes of the univariate polynomial T N + r/j^^iT 1 ^^ 1 + • • • + 770 £ k[T] 
in arbitrary order (possibly with repetitions if the polynomial is not separable). Let us 
observe that the network T evaluates the polynomials ho,. . . ,/ijv-i in the argument a 
using possibly algebraic parameters which come from the previous use of algebraic gates. 
Thus the network contains implicitly a straight -line program representation of the values 
/10(a), . . . , /iAr_i(a) which will be division free in all our applications and which contains 
possibly algebraic parameters. If T is an arithmetic network which contains algebraic 
gates we shall say that T has parameters in k . 



3 Compressing Straight— line Programs 

In this section we show how we can compress in the algorithm underlying Theorems |l|, 
§ and H| some of the straight-line programs which compute intermediate results. The 
principal outcome is the following statement : 

Lemma 1 (the afRne case) Let be given polynomials /1, •-.,/» £ k[Xi, . . . ,X n ] and 
suppose that fi, - ■ ■ , fi are represented by a division free straight-line program (3 in 
k[Xi,...,X n ] of length and depth L and £ respectively. Assume that the polynomials 
fl, . . . , fi form a regular sequence in k[Xi, . . . , X n ] and that they generate a radical ideal 
I ■= (fi, ■■■,fi)- Let W := V(fi, . . . , fi) = V(I) be the affine variety defined by /1, . . . , /» 
and denote by 5 := degW the (geometric affine) degree of W and by r := n — i its 
dimension (observe that by assumption the ideal I is unmixed and the variety W is 
equidimensional) . Suppose that Xi, . . . , X n are in Noether position with respect to the 
variety W , the variables X\, . . . , X r being free and assume that there is given by its coef- 
ficients a nonzero linear form u £ k[X r+ i, . . . ,X n ] which represents a primitive element 
for I (see Section \4-$ below). The equations fi, ■ ■ ■ , fi, the variables X\, . . . , X n and the 
linear form u determine uniquely the following mathematical objects : 
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• the minimal polynomial q £ k [X\ , X r ,u] of u modulo the ideal I . This polyno- 
mial is monic and without loss of generality separable with respect to u and satisfies 
deg u q = degq < 5 . Let us call 5 := deg u q = degq . 

• the (nonzero) discriminant p £ k[X\, . . . ,X r ] of q and polynomials iv+i, • • • ,v n 6 
k[Xi, . . . , X r ,u] satisfying the conditions 

I p = (q(u),pX r+1 - v r+1 (u), . . .,pX n - v n {u)) p 

max{deg u Vj ; r < j < n} < 5 < 5. 

Finally we assume that the polynomial p and the coefficients of q and v r +i, . . . ,v n 
with respect to u are given by a division free straight-line program /?' in k[Xi, . . . ,X r ] 
having length and depth A and A respectively. 

Under these assumptions there exists an arithmetic network with parameters in k 
which from the input circuits j3 and [3' constructs a division free straight-line program 
7 in k[X\, . . . , X r ] of length 0((i 5 + L)5 11 ) and depth 0((log2i + £)log2$) such that 7 
represents p and the coefficients of q and v r +i, . . . ,v n with respect to u. The size and 
depth of this network is 0((i 5 + L)5 11 ) + A and 0((log2i + £)log25) + A respectively. 

For the proof of Lemma [l] we need the following fact : 

Remark 1 (the affine case) Let assumptions and notations be as in Lemma 0. Then 
the degrees of p and of the coefficients of v r+ i, . . . ,v n with respect to u are bounded by 
2{degqf = 25 3 , i.e. by 25 3 . 

We show this remark later. 
Proof of Lemma [J. 

Let us introduce the following notations : 

B:=k[X 1 ,...,X n ]/L ; A:= k[X u . . . , X r ] 
K:=k(X 1 ,...,X r ) ; B' := K[X r+1 , . . . ,X n ]/(I), 



where (I) denotes the ideal generated by the set I (or by f\, . . . , /j ) in K[X r+ i, . . . , X ni 

Recall that the polynomials fi,...,f% form a regular sequence in k[Xi, . . . , X n ] and 
that the variables X\, . . . ,X n are in Noether position with respect to W , the variables 
X±, . . . ,X r being free. Taking this into account, we write tt : W — ► A r for the finite 
and surjective morphism of affine varieties induced by the coordinates X\, . . . , X r . Fur- 
thermore we deduce from our assumptions (in particular from the assumption q separable 
with respect to u) that the finite dimensional K -algebra B 1 is unramified and that 
dirnxB' = deg u q < degW = 5 holds (p5f|). From this we infer that the jacobian 



A:= det 



dfjA 

dXj J i<k<i 



r-\-l<.j <n 

is a nonzero divisor modulo /. Let us observe that the polynomial A can be evaluated 
by a division free straight-line program in k[X±, . . . , X n ] of length 0(i 5 + L) and depth 
0(log 2 i + £). 
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Let [x E k[X\, . . . , X r ] be the constant term of the characteristic polynomial of A 
modulo / . Since A represents a nonzero divisor of B we conclude \i 7^ . Furthermore we 
observe that p can be evaluated by a division free straight -line program in k[Xi, . . . , X r ] 
of length (i5L)°^ + A and depth 0(log2{iS) + £) + A , and so does the product p ■ p. 

Using a correct test sequence (see j32f| ) we are able to find in sequential time (i5L)°^ A 2 
and parallel time 0(log2(i5) + £) + A a "rational" point r] = (rji, . . . ,rj r ) E k r which sat- 
isfies (p • fx)(rj) 7^ (thus we have p{rj) / and ^(77) 7^ 0). From p[rj) / we deduce 
that the morphism it is unramified in the point r\ what implies that 1 (77) consists of 
exactly 5 nonsingular points of W . To be more precise, let 7r _1 (^) = {£1, . . . C 
be the set of these points. Then we infer from p(rj) 7^ that A(£j) 7^ holds for any 
1<1<6. 

^From 7r(^) = 77 = (rji, . . . ,rj r ) we conclude that the point has the form ^ = 
(771 , . . . , rj r , • • • , Cn' > ) £ k n where the first r coordinates 771, . . . , r/ r are rational and 

independent of the index / and the last n — r coordinates £r+u • • • > £n are algebraic (i.e. 
belong to fc and not necessarily to k ) and dependent on Z . Let the linear form u have 
the form u = 9 T+ iX r+ i + • • • + 9 n X n with 9 r+ ±, . . . , 9 n E fc . For 1 < I < 5 we write 
Q := = ^r+i^r+i + • • • + CnCn (observe that (j is typically an algebraic element of 
k). Let us consider u as an indeterminate over k and over . . . ,X r ] . In this sense 

q = q(Xi, . . . ,X r ,u) becomes a polynomial in it with coefficients in k\X\, . . . ,X r ] and 
q(rj,u) := q(r]i, . . . ,rj r ,u) becomes a polynomial in u with coefficients in k (i.e. we have 
q(ij,u) E k[u] ). Since g is monic in u and of degree 5 we conclude that q{"q : u) is monic 
and of degree 5 too. Furthermore p(r/) is the discriminant of q(rj,u) and from p(rj) we 
deduce that q(rj, u) has exactly 5 zeroes in k which are all distinct. Let us analyze this fact 
more in detail : for any polynomial p E k[X±, . . . , X n ] let us write p(rj, X r +i, . . . , X n ) := 
p(r)i,...,r) r ,X r+1 ,...,X n ) (thus we have p(r], X r+1 , . . . , X n ) E k[X r+1 , . . . , X n ] ). We 
observe that {£1 = ix~ l {rj) = {r]}xV(fi(r], X r+1 , . . .,X n ), . . . , /*(//, X r+1 , . . .,X n )) 

holds. From I p = (q(u),pX r +i — v r +i(u), . . . , pX n — v n {u)) p we deduce that for any 
1 < I < S the identities 

t, - L „ ^+lfa,Ci) Vn(T],Cl) 

and q(r], Q) = hold (as before we write Vj(rj, u) := Vj(rji, . . . ,r/ r ,u) E k[u] for r < j < n ). 
Since the points £1, • • • , £5 are all different we conclude that the values £1, . . . , Q E k are 
distinct zeroes of the polynomial q(rj,u) E k[u] . Since q(rj,u) is monic and of degree 5 
we conclude now that 

8 

q(v,u) = ]J(u - CO 
1=1 

holds. Thus the values Clj • • • > Cs are exactly the zeroes of q(ij, u) . 

Replacing in all polynomials fi, ■ ■ ■ , f%, p,q, iV+lj • • • } v n the variables Xi, . . . , X r by 
the new ones X\ — rji, . . . , X r — r} r and leaving the remaining variables X r+ i, . . . , X n and 
u unchanged, we may assume without loss of generality that r/ = (0, . . . , 0) E k r holds. 
Let 1<1<6. Recall that £, = fai, . . . , 7/r, . . . , $>) = (0, . . . , 0, £%, . . . , $>) , 

= /ite) = eS-n • • • » = /i(o, . . . , 0, eS l5 . . . , d ] ) 
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= mi) = /ifa Clt, ■ ■ ■ , = fi(0, . . . , 0, . . . , £W) 

and 

+ A(fi) = A(r/, CS 15 ■ ■ ■ , d°) = A (0> • • • , 0, • • • , d°) 

holds. Thus the point (f^, . . . , $) e fc* is a nondegenerate zero of the equation system 
given by the polynomials fi(r], X r+1 , . . . , X n ), ... , fi(r], X r+1 , X n ) £ k[X r+1 , ...,X n ] 
which are in fact the polynomials /i(0, . . . , 0, X r+ i, . . . , X n ), . . . , /j(0, . . . , 0, X r +%, . . . , X n ) . 
^From Hensel's lemma (which represents a symbolic version of the Implicit Function Theo- 
rem) we deduce that there exist formal power series RP +1 ,...,E^ £ k[[X u ...,X r }\ with 

ifl(r?) = • • • , R$(v) = ^ such that for fl« := (X u ...,X r , R^ +1 , ...,R$) the 

identities 

h(RW) = 0,...,fi(R®) = (1) 

hold in . . . , X r ]] (see |35|] or [pql). Our next task is the construction of a "sufficiently 

accurate rational approximation" of the solution of the system given by the poly- 
nomials fx, . . . ,fi. This rational approximation will be represented by rational functions 
R^ +1 , . . . , R$ € ,...,X r ) which are all defined in the point rj = (0, . . . , 0) and can 
therefore be interpreted as power series belonging to fc[[-Xi, . . . , X r ]\ . For each r < j < n 
the power series interpretation of the rational function Rj will satisfy the congruence 
relation 

if -if €(*!,..., X r ) 2 * 3+1 (2) 

where (Xi, . . . , X r ) 2p+1 denotes the ( 28 3 + 1 )-th power of the ideal generated by X\ , . . . 
... ,X r in . . . , X r \] . The i -tuple of rational functions (ifi; • • • j Rn) is defined by 

applying to the system /,..., fi at least 2>\log2&\ +2 , i.e. roughly 3|7o<?2#l +2 Newton iter- 
ation steps starting from the particular nondegenerate solution £j = (0, . . . , 0, Q+i, ■ ■ ■ , ffl) £ 
k n (here we consider /i, ■••,/« as z-variate polynomials depending on the variables 
X r+ i, . . . , X n ). In order to compute (^+i> • • • > ) we have to evaluate the polyno- 
mials ff. and their partial derivatives for 1 < /c < j and r < j < n at 3[7o<?2<5] + 2 
iteration points. By Q and |Q] (see also ||49||) this can be done by a division free straight- 
line program of size 0{LlogiS) and depth 0(tlog<i$) . Moreover we have to invert the 
jacobian matrix f^jfM i<fc<i a * the same 3|7o<72<5] + 2 iteration points. This costs addi- 

r <Cj <n 

tional 0(i 5 log25) nonscalar arithmetic operations organized in depth 0(log2ilog25) (see 
H; observe also that these matrix inversions require divisions). Thus (if_i, . . . , i?n ) are 
represented by a straight-line program ji in . . . ,X r ) which contains as intermedi- 

ate results only rational functions which are defined in the point r/ = (0, . . . , 0) . The size 
and depth of 7/ are 0((i 5 + L) log 2 S) and 0((Zog 2 2 + i) log 2 5) respectively. 

For each 1 < I < 6, let us write £W := (X u . . . , X r , ij® x , . . . , R$) . We are going to 
consider «W := «(Jj(0) = r+:L if. x + ■■■ + 9 n R^ and fi« := u(R®) = ^+1^+1 + • ■ ■ + 

9 n R^ . Observe that i^]^) = £li , . . . , 1$ (rj) = & implies u^{r]) = &. Without 
any nonscalar extra cost we may assume that 7/ computes also the rational function 
{j(0 = 9 r+1 R^ 1 + • • • + 8 n Rn^ ■ Let 7 be the straight-line program in . . . ,X r ) 



14 



which we obtain by joining all circuits 71, . . . ,75 . The straight -line program 7 computes 
. . . , and has size 0((i 5 + L)5 log 2 5) and depth 0{{log2i + £) log 2 5) . Observe 
that all the intermediate results of 7 are defined in the point rj = (0, . . . , 0) . 

Let 1 < I < 5 . Taking into account that q is the minimal polynomial of u modulo the 
ideal I = (/i> . . . , fi) we deduce from ([!]) and the assumptions of the lemma the identity 
q(X\, . . . , X r ,u^) = . Moreover (^) implies that the congruence relation 

u (0_ fi (0 e(X u ...,X r f 3+1 (3) 

holds in the power series ring . . . , X r }] . 

Recall that we consider 5 as a polynomial in u with coefficients in k[X\, . . . , X r ] . 
Moreover q is monic in u and of total degree 5 . Thus the polynomial q has the form 
q = QmU m with q m G k[X±, . . . ,X r ] and degq m < 5 for < m < 5 . Furthermore 

0<m<<5 

we have q. = 1. On the other hand we know already that q(X\, . . . ,X r ,u^ 1 ') = holds 
for any 1 < I < 5 . From (77) = Q we deduce that the power series , . . . , u'^ are 
all distinct. Since q is monic and of degree 5 in u and . . . are distinct zeroes 
of q we conclude that in • • • , X r ]][ti] the identity 

q= J] (u-uW) (4) 
1<Z<5 

holds. For < m < 5 let us denote by a m the m-th elementary symmetric function in 
5 arguments. From (0) we deduce q m = {—lY~ m a m {u <yl \ . . . ,u^) . In combination with 
(|3|) this implies that for < m < 5 the congruence relation 

q m ~ (-l) S - m a m (uV,. . . ,#)) G (X 1: . . .,X r )~ 8+1 (5) 

holds in the power series ring . . . , X r ]] . As for < m < 5 the polynomial 

q m G k[Xi, ... ,X r ] has degree at most 8 we deduce from (||) that the power series ex- 
pansion of {-lf- m a m (vP-\...,u^) G k[[Xx,...,X r ]\ coincides up to de gree 5 with 
q m . Combining the straight -line program 7 which computes the rational functions 
, . . . , with a fast and well-parallelizable algorithm for the evaluation of the set 
of elementary symmetric functions {a m ; < m < 5} (see pi], Satz 3.1. and [36], 



Chapter 8, Exercise 8.15) we obtain an arithmetic circuit 70 in k(Xi, . . . ,X r ) which for 
< m < 5 computes all rational functions (— l) s ~ m a m (u^ 1 \ . . . , u^) . The circuit 70 
has size 0((i 5 + L)Slog 2 5) and depth 0((log2i + £) log 2 S) and all its intermediate re- 
sults are rational functions of k(Xi, . . . , X r ) which are defined in r/ = (0, . . . , 0) . Taking 
this last observation into account we apply to the circuit 70 the well-parallelizable Ver- 



meidung von Divisionen technique contained in the proof of [40], Proposition 21. (With 
respect to parallelism we observe here that only one division by a suitable power of A 
is necessary at the very end of the procedure. If one is interested only in the sequential 
aspect of this technique one may apply directly the simpler algorithm underlying pCf , 
Satz 2). In this way we obtain a division free straight-line program 71 in A;[Xi, . . . ,X r ] 
which for < m < 5 computes the power series expansion in . . . ,X r ]] of all ra- 

tional functions (— l) s ~ m a m (u^ 1 \ . . . , u^ s ') up to degree 6. Since for < m < 5 these 
truncated power series expansions coincide with the polynomials q m G k[X±, . . . ,X r ] we 
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conclude that 71 represents these polynomials and hence the coefficients of q in the rep- 
resentation q = QmU m ■ The size and depth of 71 are 0((i 5 + L)5 S log 2 5) and 

0<m<<5 

0((log2i + i) log 2 5) respectively. 

Once given the straight -line program representation 71 of the coefficients q m of the 
polynomial q we are able to compute the discriminant p using additional 0(<5 5 ) arithmetic 
operations organized in depth 0(log25) . Therefore we can extend the circuit 71 to a 
division free straight-line program 72 in k[X\, . . . , X r ] of roughly the same size and 
depth as 71 , such that 72 computes also the discriminant p . 

In order to finish the proof of Lemma |l] we have to find a straight-line program rep- 
resentation of the coefficients of the polynomials v r +\, ■ • • , v n with respect to the variable 
u such that the resulting circuit has about the same size and depth as 72 . By the way 
we shall assume without loss of generality that the circuit 71 (and hence 72 ) contains the 
circuit 7 . 

Let r < j <n and let Vj = J2 OmU m with dm G k[X\, . . . , X r ] .From ([[]) and the 

0<m<<5 

assumptions of the lemma we deduce that for any 1 < I < 5 the identity 

Vj (u^)-pRf =0 (6) 

holds. Taking into account the degree bound deg u vj < 5 we interpret the identities (||) as 
an inhomogeneous 5x5 linear equation system for the coefficient vector (dg , . . . , a ^\) of 
Vj . The matrix of this equation system is the 5x5 Vandermonde matrix corresponding to 
the 5 (distinct) power series (vP 1 ' , ■ ■ ■ , u^) and its inhomogeneous part is the row vector 
we get transposing the 5-tuple (pRj, . . . , pRj) ■ From (|2|) and (||) and the assumptions 
of the lemma we deduce as before that for any 1 < I < 6 the congruence relation 

Vj {u®)-pRf e(x 1 ,...,x r f 3+1 (7) 

holds in fc[[-X"i, . . . ,X r ]] . We interpret now the congruence relations (^) as an inhomoge- 
neous 5x5 linear equation system whose (unique) solution is a rational approximation 

to the coefficient vector (oq , • • • ) a j_i) °f v j ■ More precisely we extend the circuit 72 

which computes p and for 1 < I < 5 the rational functions vP^ and Rp to a straight- 
line program 7^ in k(Xi, . . . , X r ) as follows : for 1 < I < 5 the circuit contains as 

intermediate results the rational functions pRj ■ Moreover 7^ computes rational func- 
tions 5q , . . . , which are obtained by multiplying the 5 -tuple (pRj , . . . , pRj) by 
the inverse of the Vandermonde matrix corresponding to the (distinct) rational functions 
(u^\ . . . ,u( S ') ■ We observe that all intermediate results of are rational functions of 
k(X\, . . . , X r ) which are defined in the point rj = (0, . . . , 0) . Moreover 7'- has size and 
depth 0({i 2 + L)5' 3 log 5 + 5 5 ) and 0((log 2 i + i)log 2 5) respectively. From © and © we 
deduce that for any < m < 5 the congruence relation 

aW-a«e(Ii,...,I r f +1 (8) 

holds in . . . , X r }] . 
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By Remark [Tl the degree of any polynomial ajj G k[X\, . . . , X r ] is bounded by 
2(5 3 . Therefore (||) implies that the power series expansion of the rational function a™ 
in . . . , X r ]] coincides up to degree 25 3 with a m . We join now all the circuits 

7^. 1 1, • • • > 7n an d apply to them the same Vermeidung von Divisionen technique as before. 
In this way we obtain a division-free straight-line program 7 in k[X±, . . . , X r ] which 
computes for r < j < n and < m < 5 the power series expansions in fc[[Xi, . . . , X r ]\ of 

all rational functions a m up to degree 25 . Since for r < j < n and < m < 5 this trun- 

(7) 

cated power series expansions coincide with the polynomials dm , the circuit 7 represents 
all coefficients with respect to u of the polynomials v r+ i, . . . ,v n . Without loss of general- 
ity we may assume that the straight-line program 7 extends the circuit 72 . Thus finally 
the straight-line program 7 computes p and all the coefficients of q and 77+1,- •• ,v n . 
with respect to u . The size and depth of 7 are 0((i 5 + L)5 9 log25 + 5 11 ) = 0((i 5 + L)6 ) 
and 0((log2t + i)log25) respectively. ■ 



We shall make use of Lemma |l| in the proof of Theorem [j] and Theorem ^ (i) which 
deal with the zerodimensional algorithmic elimination problem in the affine case. For the 
same problem in the toric case, i.e. for Theorem ^ and Theorem ^ (ii) , we need a slightly 
different version of Lemma |l|. 

This is the content of the next statement. 



Lemma 2 (the toric case) Let be given polynomials f\, . . . , £ k[Xi,...,X n ] and 
suppose that fi,-..,f% are represented by a division free straight-line program (3 in 

n 

k[X\, . . . ,X n ] of length and depth L and £ respectively. Let also g '■= Yl Xj . Assume 

3=1 

that the polynomials fi, ■ ■ ■ , fi form a toric complete intersection ideal I := . . . , fi) 
whose localization I g is radical. Let W* be the union of the toric irreducible compo- 
nents of the affine variety V(/i, • • • , fi) = V(I) defined by and denote by 
5* := degW* = deg*V(fi, . . . , fi) the (geometric) degree of W* and by r := n — i its 
dimension (observe that by assumption the localized ideal L g is unmixed and the "toric" 
variety W* is equidimensional) . We assume that W* is not empty, i.e. 5* > . Suppose 
that X\, . . . ,X n are in Noether position with respect to the variety W* , the variables 
X\ , . . . , X r being free and assume that there is given by its coefficients a nonzero linear 



form u S k[X r+ i, . . . ,X n ] which represents a primitive element for L g (see Section £J5 be- 
low). The equations f\,...,fi, the variables X\, . . . , X n and the linear form u determine 
uniquely the following mathematical objects : 

• a polynomial q* £ k[X\, . . . ,X r ,u] which is monic and without loss of generality 
separable with respect to u and satisfies deg u q* = degq* < 5* 

• the (nonzero) discriminant p* £ k[X\, . . . , X r ] of q* and polynomials . . . 6 
k [X\ , X r , u] satisfying the conditions 

Ip*9 = (q*( u )iP* X r+l ~ W* +1 (V), . . .,p*X n - V* n (u)) p *g 

max{deg u v* ; r < j < n} < degq* < 5* . 
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Finally we assume that the polynomial p* and the coefficients of q* and v*<i, . . . , 
with respect to u are given by a division free straight-line program (3* in k[X\, . . . ,X r ] 
having length and depth A* and A* respectively. 

Under these assumptions there exists an arithmetic network with parameters in k 
which from the input circuits (3 and [3* constructs a division free straight-line program 7* 
in k[Xi,...,X r ] of length 0((i 5 + L)5* u ) and depth 0((log 2 i + i)log 2 S*) such that 7* 
represents p* and the coefficients of q* and v* +1 , . . . , v * with respect to u . The size and 
depth of this network is 0((i 5 + L)5* n ) + A* and 0((log2t + i)log20~*) + A* respectively. 
This statement remains true if "toricity" is understood with respect to a linear change 

n 

of the original variables X% , • • • , X n and the polynomial g = Yl Xj is replaced by the 
product of the new variables. 

Since the proof of Lemma [l] is completely local, only minor changes at its very beginning 
and the application of a suitable toric version of Remark [l] at its very end are necessary 
to convert it into a proof of Lemma ||. The essential steps and arguments are almost 
textually the same. The most subtle point is the toric version of Remark || we need at the 
end of the proof. For the statement of this toric version we replace just in Remark |l] the 
polynomials p , v r+ ±, . . . , v n by p* , v* +1 , . . . , u* and the degree bound by 2(degq*) 3 , i.e. 
by 2(5* ) 3 . We shall come back to the proof of that at the end of this section. In order to 
avoid repetitive argumentation we omit the proof of Lemma ^. 

The following Remark |^ refers to work in progress. It answers the question whether 
computing with algebraic elements of k is really necessary in the algorithm underlying 
Lemma |l| and ||. On the other hand one might ask whether to this algorithm corresponds 
a counterpart in the bit model of boolean circuits in case k:= Q . We postpone the answer 
to this second question to the end of the paper (see Section [?]). 

Suppose now that k is a hilbertian field for which a factorization algorithm for uni- 
variate polynomials is available at "moderate arithmetic costs". This means that we may 
add (in a computationally reasonable way) to our arithmetic networks over k special fac- 
torization gates for univariate polynomials of degree D , where D is an arbitrary (but 
fixed) natural number. Our new complexity model takes such a gate into account at costs 
of with respect to sequential time (this is quite realistic in view of |4j|) and at 

costs of 0(log2D) with respect to parallel time (this is rather cheap). We call such a 
base field k hilbertian with univariate polynomial factorization at moderate costs and the 
corresponding algebraic complexity model is called arithmetic network ( over k ) with fac- 
torization gates. Observe that this model is nonuniform by the way we have introduced it. 
However in view of |2£| and the evidence that "efficient" (rather than "effective" ) versions 
of Hilbert's Irreducibility Theorem must exist (see e.g. [l4j| , [pcfl , f5S( l and |6^| ), we may 
hope to obtain at least a reasonable randomization (if not uniformization) of the algorithm 
underlying the next statement which we pronounce without proof. 

Remark 2 Suppose that k is a hilbertian field with univariate factorization at moderate 
costs. Suppose that there is given the same situation (with the same notations) as in 
Lemma [1] and || and make the additional assumption that the circuits (3' and (3* use only 
parameters from the base field k . Then, subject to the following modifications, the same 
conclusions as in Lemma and are true : 
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• the output network in question is an arithmetic network over k with factorization 
gates 

• the straight-line programs 7 and 7* are division free circuits in k[Xi, . . . , X r ] 
(which use only parameters from k ). 

We finish this section showing the assertion of Remark [l] and of its toric counterpart. 

Proof of Remark [J. Let assumptions and notations be as in Lemma |l| and its proof. 
We shall consider the A: -algebras C := k[X\, . . . , X r , u]/q(u) and 

C p :=k[X l ,...,X r ] p [u\/{q{u)) p 
and compare them to the k -algebra B := k[X\, . . . , X n ]/(fi, . . . , fi) and its localization 

Bp := k[X\, . . . , X r ] p [X r+ i, . . . , X n )/ (fi, . . . , fi) p . 

Let us denote by ~ the canonical k -algebra homomorphism which associates to each ele- 
ment in . . . , X n ] , 

k[Xi, . . . , X r ] p [X r+ i, . . . , X n ] , k[X\, . . . , X r , u] or k[Xi, . . . , X r ] p [u] its residue class in 
B , Bp, C or C p . Observe that the condition (/1, ■ ■ • , /j)p = I p = (q(u),pX r+ i — 
v r+ i(u), . . . , pX n — v n {u))p induces a A; -algebra isomorphism r : B p — > C p which leaves 

k[Xi, . . . ,X T ]p fixed and which satisfies the condition T(X r+ \) = r+1 ^ , . . . ,r(X n ) = 

P 

— . Since the polynomials f\,...,f% form a regular sequence in k[Xi, . . . , X n ] and the 

P 

variables X±, . . . , X n are in Noether position with respect to the ideal (fi, ■ ■ ■ , f%) , the 
A; -algebra B is a free k[X±, . . . ,X r ] -module of rank 5 (see e.g. [pq| , Lemma 3.3.1). The 
same is true for the A; -algebra C . Hence B p and C p are free k\X\, . . . , X r ] p -modules 
of rank 5. The polynomial is q by assumption separable with respect to the variable 
u. Thus we conclude that B p and C p are unramified A;[Xi, . . . , X r ] p -algebras (this fact 
has already been used in the proof of Lemma |l|). Therefore the (ordinary) traces which 
map B and C onto k[X\, . . . , X r ] are nonzero and the corresponding traces of B p and 
Cp are nondegenerate (observe that p is the discriminant of the k[X\, . . . ,X r ] -module 
basis of C given by the elements 1, u, . . . , , u <5 ~ 1 ). Since the k -algebra isomorphism r is 
k[X\, . . . ,X r ] p -linear it leaves the traces fixed. We denote them therefore by the same 
symbol "Tr" (see [42|, Appendix F as reference for traces in noetherian algebras over 
rings). Let r < j < n . As in the proof of Lemma |l| let us write Vj = amU m 

0<rn<8 

with Om £ k[Xi, . . . ,X r ] . Observe that pXj = Vj(u) = J2 amU m holds in B. For 

0<m<5 

< m' < 5 consider the polynomial 



Si) 



pTr(X J u m ')= J2 ^Tr(u m+m ') (9) 
0<m<<5 



which belongs to k[X\, . . . , X r ] . 

Observe that Tr(XjU m ) is also a polynomial of k[X±, . . . , X r ] which satisfies by [52], 
Theorem 13 the degree bound degTr(XjU m ') < degWdeg{X jU m ') = {m' + 1)5 < P . Thus 
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S3) 

cfj is divisible in k[X\, . . . ,X r ] by p and satisfies the degree bound 



P 
Jj) 

deg^- < 5 2 . (10) 

In the same way we see that for any < m < 5 the trace Tr(u m+m ') is an element 
of k[X\, . . . , X r ] which satisfies the degree condition 

degTr(u m+m ') < 25(5- 1). (11) 

For j fixed and m verifying < m < 5 we interpret the equations (||) as an inhomo- 
geneous 5x5 linear equation system for the coefficient vector (o,q\ . . . , ■ ^he m at rrx 

of this equation system is [Tr(u m+m ) ) _ whose entries belong to k[X±, . . . , X r ] . 

V /0<m,m'<8 

The determinant of the equation system is therefore p, which is the discriminant of the 
k[Xi, . . . , X r ] -module basis 1, u, . . . , -u" 5 " 1 of C . The inhomogeneous vector of the equa- 
tion system is obtained transposing the <5 -tuple (cq \ . . . , c^_ A ) ■ 

^From Cramer's rule and the degree bounds © and (|ll]) we conclude therefore 



dega$ < 26 



for any r < j < n and < m < 5. Similarly we deduce from (11) the degree bound 
degp < 25 3 < 25 3 . m 



The toric version of Remark [T] can be proved essentially by the same method. One has 
only to replace the -algebra B by the (reduced) coordinate ring of the "toric" variety 
W* . Of course this coordinate ring is not a free k[X\, . . . , X r ] -module anymore, but it is 
still a finite faithful module over k[X\, . . . , X r ] , i.e. an integral extension of k[X±, . . . , X r ] . 
This is sufficient for traces to be well defined and for |3^], Theorem 13 to be true also in 
this new context. The rest of the proof of the toric version of Remark || is textually the 
same as in the affine case. 

Let us also observe that Remark [l] can be formulated with a slightly better degree 
bound of 5 2 + 1 , i.e. of 5 2 + lin the affine case and (degq*) 2 + 1 , i.e of 5* 2 + 1 in the 
toric case. This can be easily deduced from our elementary arguments in Section [5] taking 
into account that all the polynomials involved are homogeneous. 



4 Technical Lemmas 

4.1 Squarefree Representation and Greatest Common Divisor Compu- 
tation of Univariate Polynomials 

Let R be an integral and factorial k -algebra with fraction field K , T a new variable, D 
a fixed natural number and P,Q£ R[T] two polynomials of formal degree D . We think 
that P and Q are given by their coefficient vectors each of length D . In the statements 
which follow now we refer to (ordinary) arithmetic networks and straight -line programs 
over R which receive the coefficient vectors of P and Q as inputs. Furthermore greatest 
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common divisors (gcd's) will always be taken with respect to the principal ideal domain 
K[T] although they are represented in R[T] (thus they are not unique up to units in R). 
The proofs of the following well-known lemmas can be found in [40]. 



Lemma 3 There exists a division free well parallelizable arithmetic network V of size 
0(D e ) which from the coefficients of P and Q computes the coefficients of a greatest 
common divisor (belonging to R[T] ) of these polynomials. 

The greatest common divisor computed by the algorithm underlying this Lemma is a 
(not necessarily primitive) polynomial of R[T] which we denote by gcd(P, Q) . A similar 
well parallelizable ( 0(s 5 D e ) sequential time ) complexity bound holds for the computation 
of the simultaneous greatest common divisor of more than two (say s ) polynomials of 
degree at most D belonging to R[T] . 

Sometimes the second polynomial Q is given by a division free straight-line program 
in R[T] of length L and depth t instead by its coefficients. In this case we have 

Lemma 4 There exists a division free arithmetic network of size LD°W and depth 
0(log2(D) + £) which computes the coefficients of the greatest common divisor of P and 

Q. " 

The algorithm underlying Lemma ^| uses linear algebra subroutines which deal with 
square matrices over R of size at most 2D — 1 . 

As before this lemma can be generalized to the case of the computation of the simulta- 
neous greatest common divisor of more than two (say s ) polynomials, one of them given by 
its coefficient vector and having degree at most D and the others given by a division free 
straight -line program in R[T] of length L and depth I. The outcome is then a division 
free arithmetic network of size (L + s)D°^ and depth 0(log2(D) + €) which computes 
the coefficients of the greatest common divisor in question. Again the linear algebra sub- 
routines deal only with square matrices of size at most 2D — 1 . Lemma |3] can be used 
in order to compute a separable (and hence in K[T] squarefree) polynomial P £ R[T] 
which has the same zeroes as P (we call such a polynomial P a separable representation 
of P). If the characteristic of k (and hence of R and K) is zero this is a immediate 
consequence of the next lemma putting P := P* . In case of positive characteristic we 
need a more refined analysis of the situation. 

Observe that the coefficients of the derivative P' of P are immediately obtained from 
the (given) coefficients of P . Therefore we have 

Lemma 5 Suppose that the derivative P' is nonzero. Then there exists a division free well 
parallelizable arithmetic network of size 0(D e ) which from the coefficients of P computes 
the coefficients of a univariate polynomial P* £ R[T] and a non-zero element of the ring 
6 £ R, such that : 

P 

p* _ Q 

gcd(P,P') 

holds. 
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Suppose now that the characteristic of k , K and R is positive, say p . In all our 
applications of the technical lemmas of this section R will be a polynomial ring over 
k generated by some k -linear forms in the variables X\ , . . . , X n , say Y\ , . . . , Y r with 
< r < n (see Section ||). As one easily verifies there is no essential change to the 
outcome of our main algorithm in Section [B| if we replace the variables X\ , . . . , X n of our 
input equation system /i, . . . , /„ by suitable p k -th powers of them with k E IN not too 
big (recall that we suppose the ground field to be perfect and our arithmetic networks 
and straight -line programs to include special gates for the extraction of p-th. roots in 
k and in k ) . Thus we may suppose without loss of generality that we are able to extract 
suitable p -th roots in R (and hence in K ) of not too high degree (p k < D will suffice) . 
Applying now Lemma [| iteratively (with at most log p D iterations) we see that we obtain 
in sequential time 0(D 7 ) a separable representation P E R[T] of the polynomial P . 
Taking into account that in all our applications the factor appearing in Lemma || will 
be monic with respect to one of its variables and using standard tricks of parallelization 
we get a division free and well parallelizable arithmetic network of size which from 

the coefficients of P computes the coefficients of a separable representation P E R[T] of 



P (see [24], 2.1. and [26] for more details about how to modify algorithms in order to 



achieve the assumptions which allows to compute P). 

4.2 Correct Test Sequences and "Vermeidung von Divisionen" 

The algorithms in g9|, (2§|, |§, g0|, Q rely heavily on the use of "correct test sequences" 
(fl3||, [32], |40|1 ) as a tool for deciding identity of polynomials given by straight-line pro- 
grams. Unlike the probabilistic identity tests [15], p4§ , [55], [64] the choice of a suitable 
correct test sequence does not depend on the specific polynomials whose identity has to 
be checked but only on the number of variables and the size of the input circuit. The 
outcome are nonuniform deterministic or random algorithms of a stronger type as those 
proposed in [55] and [64]. 

Definition 1 Let J 7 be a set of polynomials of k[X\, . . . ,X n ] such that belongs to T . 
Let Q be a subset of k n . Q is called a correct test sequence (or questor set) for T if for 
any P E T the following implication holds : 

P(x) = for all x E Q ^ P = . 

The cardinality #Q is called the length of the test sequence Q. 

The existence of short correct test sequences is warranted by the following fact : 

Lemma 6 Let T be the class of all polynomials of k[X\, . . . , X n ] which can be evaluated 
by a straight-line program in k(X\, . . . , X n ) of (nonscalar) size L and depth i. Let 
uj ■= (2 e+1 - 2) (2 e + l) 2 and a := 6 {IL) 2 . Then for any collection O of uj 
elements of k the set Vt n C k n contains at least iv na (1 — u>~6) correct test sequences of 
length a for the class T . 

If the characteristic of k is zero we can always make the standard choice f2 := 
{1, . . . , w} C 2Z . 

Lemma ^ implies that for L > and I > there always exist correct test sequences for 
J- having length a which is polynomial in L , and that any random choice of a elements 
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a "I 

in O n leads to a correct test sequence for T with an error probability of uj~ e << ± . 
The proof of Lemma |6| can be found in |Q (see also []32|). 

We shall make frequent use of a method called "Vermeidung von Divisionen" which is 
due to and which allows to transform any straight -line program V (with parameters in 
k or k ) which computes a polynomial and contains essential divisions into an equivalent 
division free one, say V . If the size and the depth of T are L and £ respectively, then the 
size and depth of V are D 2 L and 0(£) (see [^0| for the proof of this precise statement). 



4.3 Construction of a Primitive Element 

Primitive element constructions (or "shape lemmas") are crucial for solving zerodimen- 
sional polynomial systems symbolically (see Q, @, 0], ||, |J], |3, 0, @, @, 



They furnish a handy description of the multiplication tensor of the finite dimen- 
sional k — algebras which arise as (typically reduced) coordinate rings of such equation 
systems. 

Let R be an integral k -algebra with fraction field K . Let 7 be a zero-dimensional 
ideal of K[X\, . . . , X n ] . We denote the radical of I by \fl . 

For any maximal ideal Q associated to I , let K(Q) be the field 

K(Q):=K[X 1 ,...,X n ]/Q. 

A linear form u = t\X\ + • • • + t n X n E R[X±, . . . , X n ] is said to be a primitive element 
for / (or for the finite dimensional K — algebra K[X±, . . . ,X„]/\/7) if it satisfies the 
following two conditions : 

• for any maximal ideal Q associated to I , the minimal equation (minimal polyno- 
mial) of u modulo Q has degree [K{Q) : K] sep (here [K(Q) : K] sep denotes the 
separability degree of the finite field extension K — > K(Q)). 

• For any two different maximal ideals Q and Q! associated to I , the minimal poly- 
nomials of u modulo Q and Q! are distinct. 

If the algebra K[X\, . . . , X n ]/^fl is unramified over K then a linear form u G 
R[X\, . . . , X n ] represents a primitive element for / if and only if the set {1, u, u 2 ,u 3 , . . .} 
generates modulo \fl the K— vector space K[X\, . . . ,X n ]/\/l . We shall also apply the 
notion of "primitive element" to zero dimensional localizations (by a given polynomial) of 
generally nonzero dimensional ideals. When doing so (e.g. in Lemma ||) we have to think 
/ as the intersection of all primary components of the given ideal which do not contain 
the polynomial. 

For the construction of a primitive element from a given set of generators of / we 
follow the lines of [24] or [fl(]].In order to explain our method we introduce new variables 



T% , . . . , T n and for 1 < j < n we consider the following domains : 

Rj := R[Ti , . . . , Tj_i , Tj + i , . . . , T n ] 
Kj := K(Ti, . . . , Tj_i,Tj +1 , . . . , T n ), 
and the following linear form in X±, . . . ,Xj—i, -Xy+i, ■ ■ ■ , X n : 
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Zj :— T\X\ + • • • + Tj-iXj^i + Tj + \Xj + \ + • • • + T n X n . 
Furthermore we consider the following linear form U in X\ , . . . , X n : 

U:=T 1 X 1 + --- + T n X n . 

We observe that for any 1 < j < n the identity U = Zj ■ + TjXj holds. Let T be an 
additional variable and D a natural number. The way how we shall generate a primitive 
element is contained in the following statement : 

Lemma 7 Let be given the following inputs : 

• monic polynomials fx, . . . , f n G R[T] of degree at most D such that for each 1 < 
j <n the specialization fj(Xj) belongs to the ideal I . 

• for each 1 < j < n a polynomial gj G Rj [T] which is monic in T and which has total 
degree at most D such that the specialization gj(Zj) belongs to the ideal Kj§Z) K I . 

We suppose further that fx, ■ ■ ■ ,f n and gx,- ■ ■ ,g n « r e represented by their coefficients 
with respect to the variable T. Moreover we assume that the coefficients of g\,...,g n 
themselves are given by a division free straight-line program in R[T\, . . . ,T n ] of length L 
and depth I . 

Then there exists a division free arithmetic network in R of size 

(nDL)°^ and depth 

0(log2(nD) + £) which computes the coefficient representation of the following items : 

i) a primitive element u = \\Xi + • • • + X n X n G k[Xx , • • • , X n ] for the ideal I , 

ii) a non-zero element p G R and for each 1 < j < n a polynomial Vj G R[T] such 
that pXj — Vj(u) G vJ holds, 

Hi) a monic polynomial q G R[T] such that q(u) G vJ holds. 

Moreover we have degvi < degq for 1 < % < n . From the observations at the end 
of Section J^A_ we deduce that we may suppose without loss of generality that q is sepa- 
rable with respect to the variable T and hence squarefree. This implies that the algebra 
K[Xx, • • • , X n ]/y/7 is unramified over K and that \fl = (q(u), pX\ — vx(u), . . . , pX n — 
v n (u)) holds in K[X\, . . . , X n ] . 

Lemma |^ is contained in pQ| , Proposition 27 and its proof (see also 24], Sections 3.4.6 
and 3.4.7). 



5 The Algorithm 

In this Section we describe the algorithmic procedure underlying Theorem |l], ^ and |3[ 

Let us fix the following notions and notations : let Xq be a new variable. For any 
non-zero polynomial p G k[Xx, ■ ■ ■ , X n ] of degree D we define its homogenization ^p 
(with respect to the variable Xq ) as 

h p:=X D p(Xx/X^...,X n /X ). 
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For p = we put h p:= 0. Of course h p is a homogeneous polynomial in k[Xo, . . . ,X n ] 
of degree D . For any ideal I of k[X±, . . . , X n ] , we define its homogenization h I as the 
(homogeneous) ideal in k[X , . . . , X n ] generated by the set of polynomials { h f :/£/}. 

We maintain assumptions and notations of Section [l]. In particular we suppose to 
be given input polynomials fi,---,f n € k[Xi, . . . , X n ] of degree at most d which are 
encoded by a division free straight-line program of size L and depth I . We suppose that 
fl, . . . , f n form an affine or toric complete intersection, according to the problem we are 
considering. Let H be a non-zero linear form of k\X\, . . . ,X n ] given by its coefficients. 
The fundamental problem we want to solve is the following : find a non-zero polynomial 
p (in the affine case) or a non-zero polynomial p* (in the toric case) belonging to k[T] 
such that p(H) vanishes on V = V(/i, • • • , fn) (in the affine case) or p*(H) vanishes on 

n 

V* = V(/i, . . . , f n ) \ V{ n Xi) (in the toric case). This is the content of Theorem [TJ and 

r-, i=l 

I 

^From [|^] Lemma 13 we deduce that the homogeneous polynomials h f\ , . . . , h f n can 
be evaluated by a division free straight-line program in k[Xo, . . . , X n \ of (nonscalar) size 
and depth d(d + 1) 2 L and log2d + 21 respectively. 

For 1 < i < n we denote by 1%:= (fi, ■ ■ ■ , fi) the ideal generated by /i, . . . , /, in 
k[X±, . . . , X n ] . Let us now introduce the specific notions and notations which we need in 
the affine case. Fix 1 < i < n and let Jj := \f^Ti be the radical of the homogenization 
of U = (/i, . . . , fi) . Thus Ji is a homogeneous and unmixed radical ideal of codimension 
i and the corresponding projective variety does not contain any irreducible component at 
infinity. Furthermore the homogeneous polynomial /j+i is a non-zero divisor modulo Jj 
for i < n. Let 2j:= ( h fi, Jj_i) be the homogeneous ideal generated by h fi and Jj_i . We 
observe that Jj is the intersection of all codimension i prime homogeneous ideals which 
contain Xj but not the form Xq (of course, these are associated primes of Xi). 

Let us finally consider the specific notions and notations of the toric case. Let 1 < 
i <n. We denote by J* the intersection of all homogeneous prime ideals of codimension 

n 

i which contain "Jj but not the form FJ Xi (these prime ideals are again associated to 
h Ii). By hypothesis (the family fi,...,f n forms a toric complete intersection) the variety 

n 

V(Ii) is not contained in 1^(11^) an d therefore J* ^ (Xq, . . . ,X n ) . Furthermore the 

i=l 

homogeneous polynomial is n °t a zero divisor modulo Jf t for i < n. Finally, as 

in the affine case, J 4 * is the intersection of all codimension i homogeneous (associated) 

n 

prime ideals which contain ( h fi, J*„ 1 ) but not the form FJ Xi . 

With these notions and notations fixed we are now able to describe the principal items 
which are produced as intermediate results by the algorithm underlying Theorems |l|, ^ 
and |3| and to indicate their main properties. 

For the affine as well as for the toric case the algorithm proceeds in an analogous 
manner in n recursive steps. We give now a simultaneous account of the i-step of the 
algorithm both for the affine and the toric case (here 1 < i < n). The algorithm produces 
first a linear change of the variables Xq, . . . , X n into new variables Yq, . . . ,Y n such that 
the following canonical ring homomorphism becomes a generically unramified integral 
extension : 
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in the affine case : 

k[Y ,...,Y n _ i ]^k[Y ,...,Y n ]/J i 



• in the toric case : 



k[Y ,...,Y n ^] 



k[Y ,...,Y n ]/J* 



Next the algorithm generates a nonzero linear form Ui 6 k[Y n -i+x, . . . ,Y n ] and a 
homogeneous polynomial of k[Yo, . . . , Y n _i, T] which is monic and separable with re- 
spect to the variable T and which we denote by qi in the affine and by q* in the 
toric case. The linear form m is a primitive element for the ideal Ji in the affine 
and for J* in the toric case having minimal polynomial qi or q* respectively. Finally 
the algorithm produces in the affine case homogeneous polynomials pi £ k[Yo, . . . , Y n _^\ 
and . . . , vli £ k[Yo, . . . ,Y n ~i,T] with pi ^ such that in the localized ring 

k[Y , ... , Y n -i] Pi [Tn-i+i . . . , Y n ] the ideal generated by 



is identical with the ideal (Ji) p ■ The polynomial pi will be the discriminant of the 
(separable) polynomial qi with respect to the variable T . 

A similar statement is valid for the toric case if we replace qi by q* and Jj by J* . In 
order to avoid accumulation of diacritic symbols in the proof we shall not use extra notation 

(i) (i) 

for the polynomials which play in the toric case the role of pi and v n _ i+1) . . . ,v n . The 
same notation will be applied in the affine case as well as in the toric case, no ambiguity 
will arise from that. Thus for example pi will be interpreted as the discriminant of qi or 
as the discriminant of q* , following the context. 

(i) (i) 

Let us also remark that the degree of v n _ i+1) . . . ,Vn with respect to the variable T 
will be strictly less than the degree of qi or q* respectively (this latter condition makes 

«n-i+D • • • ' v n unique in k(Y , . . . , Y n ^i)[T] ). 

Thus the k(Yo, . . . , Y n -i)— algebra k(Yo, . . . , l^_j)[T]/ (qi) (or equivalently 
k(Y ,. . . ,Y n _i)[ui]/(qi(ui))) is isomorphic to h(Y , . . . ,Y n _i)[Y n - i+ i, . . . ,Y n }/(Ji) in the 
affine case (and the isomorphism can be effectively computed). Both algebras are unram- 
ified. An analogous statement is true in the toric case. 

The linear form Ui will be given by its coefficients with respect to the variables 
Yn-i+li ■ ■ ■ jY n , whereas qi , q* , v^_ i+l , . . . , Vn will be given by their coefficients (which 
belong to k[Yo, . . . , Y n -i] ) with respect to the variable T . These coefficients and the poly- 
nomial pi will be represented by a division free straight-line program in k[Yo, . . . , Y n ^i] ). 
Furthermore our procedure will produce this straight -line program. 

Finally let us observe that these items which appear as intermediate results of our 
procedure are canonical and intrinsic objects with a precise geometric meaning (on just this 
meaning is based the data compression contained in Lemma |l] and ^ which is fundamental 
for our approach) . Let us consider only the affine case : to the linear change of the variables 
Xq, . . . ,X n into Yq, . . . ,Y n corresponds a finite surjective morphism of affine varieties 
^i ■ V(fi, . . . , fi) — ► A n_i . The linear form Uj complements 7Tj to a finite morphism 
(7Tj, Ui) which maps V(/i, . . . , fi) onto a hypersurface of A n_l+1 . The minimal equation 
of this hypersurface is just the polynomial we obtain specializing in qi £ k[Xo, . . . ,X n ] 
the variable Xq to 1 . With the same type of specialization we see that pi describes 



qi(ui),piY r 




n— i+1 



(u^, . . .,piY n - v$(v,i) 



26 



the discriminant variety of the generic fiber of 7Tj and that the polynomials f„-i> ■ ■ ■ i v n 
represent a "univariate" rational parameterization of the ni -fibers of the points of A n_i 
which lie outside of this discriminant variety. 

5.1 The Recursion 

Keep 1 < i < n and notions and notations of the last subsection fixed. We are going 
to describe the i + 1-th step of our recursive main algorithm. Let Si := degV{f\, ...,/») 
and 5* := deg*V(fi, . . . , fi) and recall that degfj < d holds for 1 < j < n. With 
these notations we have degxqi = degqi < 5i , degxq* = degq* < 5* . In order to simplify 
notations we shall assume from now on that the polynomials qi and q* have degree exactly 
5{ and 5* respectively. This assumption does not restrict the generality of our arguments 
(compare Section Let Tj be a division free straight -line program in k\Yo, . . . , Y n ~i] 
which represents the coefficients of Uj and which evaluates pi and the coefficients with 
respect to T of the polynomials v^}_ i+1 , . . . ,Vn and of qi in the affine, and of q* in the 
toric case. Denote by Aj and Aj the size and depth of respectively. 
The core of our procedure is the following technical result : 

Proposition 1 There exists an arithmetic network with parameters in k of 

• size i(d5iLki) ^ 

• depth 0(log 2 (d5i) + £) + A; 

which from T{ as input produces a linear change of variables Yq, . . . ,Y n into new variables, 
say Yq, . . . , Y^ , such that Yq, . . . , Y£ are in Noether position with respect to the ideal Jj+i 
in the affine case and to J* +1 in the toric case, the variables Yq, . . . , Y^ l _ i _ 1 being free. 

Furthermore the network produces the coefficients of a linear form in the ring 

M^n-ii • • • ' ^n] an d a division free straight-line program in kY^, . . . , Y^^i] which rep- 
resents the polynomial pi+\ and the coefficients with respect to T of , . . . ,Vn^ G 
k[Y \ l^_ i _ 1 ,T] and of q i+1 or q* +1 G k[Y ', Y^_ i _ 1 ,T] , following the (affine or 
toric) case. 

The size and depth of this straight-line program are i(L+Ai)(d5i)°^ and 0(log2(d5i)+ 
£) + Aj respectively. The parameters of the arithmetic network and of the straight-line pro- 
gram it produces are contained in the field extension of k generated by the parameters of 
Ti. 

We observe without proof that the algorithm underlying this procedure can be or- 
ganized in such a way that it uses only linear algebra subroutines dealing with square 
matrices of size at most 2d5i (or 2d5* ) . 

We divide the i + 1 -th recursive step of our algorithm in three parts : 

• recursive Noether normalization. 

• recursive generation of a primitive element. 

• cleaning extraneous irreducible components. 
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The first two parts are common for both the affine and the toric case. Only the third 
part distinguishes in some few technical points between the two cases. 

Before giving the details of these three parts, let us recall the following procedure to 
compute resultants which is implicitly contained in P| and [52| : 



Lemma 8 Let J be an unmixed homogeneous radical ideal of k[Xo, . . . ,X n ] and let us 
suppose that the canonical homomorphism 

<t> : k[X , X n ^} — k[X , . . . , X n ]/J. 

is an integral graded ring extension. Then, for an homogeneous polynomial g in the ring 
k[Xo, . . . , X n ] its minimal polynomial m g £ k[Xo, . . . , X n -i, T] has the form 

T D + a D -iT D - 1 + ■ ■ ■ + a , 

where D < degJ and aj £ k[Xo, . . . , X n _i] is zero or a homogeneous polynomial of degree 
(D - j) ■ deg(g) . 

The same conclusion remains true if we replace the minimal polynomial m g by the 
characteristic polynomial Xg °f 9 ■ Moreover the minimal polynomial m g is squarefree if 
4> is generically unramified and we have D < degV(J) in this case. 

Proof. Observe first that our assumptions on J imply that m g and Xg belong to 
k[X ,...,X n „i,T]. 

N 

Let J = P| Qj , where the Qj 's are the homogeneous prime ideals of codimension 

3=1 

i associated to J. From our assumptions we obtain for 1 < j < N an integral ring 
extension : 

fa : k[X , X n _i] ^ k[X Q , X n ]/Qj 

which is generically unramified if the same holds for (j) ■ We consider the multiplication by 
the polynomial g modulo the ideal Qj as a k[Xo, . . . , -linear map. In this sense we 
shall speak in future about "the homothety defined by multiplication by g modulo Qj " . 
The minimal polynomial m g of this homothety verifies the conclusions of the lemma with 
respect to degrees (see ]|], Remark 9 and its proof). Moreover fa generically unramified 
means that the corresponding finite field extension is separable. Therefore the minimal 
polynomial m g is squarefree if eft is generically unramified. Since the minimal polynomial 
of the homothety defined by the multiplication by g modulo J is a product of some of the 
given rrig , • • • , rrig (without repetitions if (j) 1S generically unramified), the first assertion 
of the lemma follows. The third one is a consequence of the fact that the polynomials 
vfig 1 ^ , . . . , Trig are irreducible and that J is radical. The second assertion follows from 
the observation that the irreducible factors of Xg are the polynomials m g , . . . , rrig . ■ 

Note that if g is not a zero divisor modulo J , the term oo is a non-zero homogeneous 
polynomial of k[Xo, . . . , X n _j] of degree D ■ deg(g) . 
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5.1.1 Recursive Noether Normalization. 

Without loss of generality we describe just the affine case (the toric case follows simply 
by replacing Jj by J* and so on). 

First we observe that fi+i is not a zero divisor modulo Ji . According to the notations 
of Lemma |8| above let r/uj be the minimal polynomial of the homothety given modulo 
Ji by h f i+ i : 

m hfi+1 = T D + gd^T - 1 + • • • + a . 

The coefficient ao £ k[Yo, . . . , Y n -j\ is a non-zero homogeneous polynomial of degree at 
most d5i and x is squarefree by Lemma |8| 

In order to obtain a Noether normalization with respect to = ( fi+i, Ji) (and 
hence with respect to Jj+i ) we have just to find a k -linear change of the variables 
Yq, . . . , Y n -i into new ones, say Y \ . . . , Y! n _ i , such that the (uniquely defined) polynomial 
A (Yq, • • • > ^n-i) e k[Y^, . . . , = k[Y , . . . , Y n ^i] which satisfies the equality 

A (Y^...X-i) = ao(Y ,...,Y n „ i ) 

is monic in the variable Y^_j • We obtain this linear change of variables as follows : 
First observe that modulo Ji the forms 1, Uj, . . . , uf-' 1 represent a basis of the 
k(Y , ... , Y n _i) -vector space k(Y , ... , Y n _i)[Y n _ i+1 , ... , Y n ]/(Ji) . By assumption the 

(i) (i) 

coefficients with respect to the variable T of the polynomials qi , v n _ i+1 , . . . , v„. £ 
k[Yo, . . . , Y n ^i, T] as well as the polynomial pi £ k[Yo, . . . , Y n ] are given by the circuit 
Ti . This implies that we can write down directly the (homogeneous) entries of the com- 
panion matrix Mj £ k[Yo, . . . , Y n _i] SiXSi of qi . The matrix Mi describes the homothety 
defined by the multiplication by Ui in k(Yo, . . . , Y n _i)[Y n ^i + i, . . . , Y n ]/( Ji) with respect 
to the basis given by the forms 1, Ui, . . . , uf 1 ^ 1 . With respect to the same basis the ho- 
mothety defined by the multiplication by h fi+i has a matrix representation Mhf i+1 which 
satisfies the equation : 

tfM hfi+1 = h f l+1 (Y , . . . , Y n . u vf_ i+1 {Mi), . . . , t,« (Mi)) 

for some k < d (in a slight abuse of notation we write h fi+i{YQ, • • • ^Y n ) for the polyno- 
mial obtained from the original h fi + i(Xo, . . . ,X n ) by means of the change of variables 
from Xq, . . . , X n to Yq, . . . ,Y n and multiplying in this new expression submonomials in 
Yq,... , Y n ^i by suitable powers of pi ). Observe that p\Mh^ i+l £ k[Y , . . . , Y n _i] 5iXSi holds 
and that the entries are homogeneous polynomials. Therefore the entries of the matrix 
pfMhf i+1 can be computed by an (ordinary) straight-line program in k[Yo,... ,Y n ^i] of 

length Ldf^ and depth 0(log2di) from the coefficients of qi , v^_ i+l , . . . , Vn and 
from pi . This is simply done by applying the evaluation scheme for h fi+\ (which is given 
by the input) to suitable entries of the matrices v^_ i+l (Mi), . . . ,Vn\Mi) which replace 
the variables Y n _ j , . . . , Y n . 

Next, we use Berkowitz's well parallelizable and division free algorithm || in order 
to compute the (homogeneous) coefficients of the characteristic polynomial X p khj i+1 £ 

k[Yo, . . . , Y n _i, T] of the homothety given by p^ h fi + i modulo Ji . The computation of the 
coefficients of X p k hf i+1 from the entries of p\Mh^ i+l requires an additional division free 

and well parallelizable straight -line program in k[YQ, . . . , Y n _i\ of size . 
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Now taking into account our observations at the end of Section 4.1 we apply Lemma 
H in order to obtain the (homogeneous) coefficients with respect to T of a separable 
representation x kh f^ ^ k\Yo, • • • > Y n -i, T] of x «f , , • This separable representation 
is without loss of generality the minimal polynomial rn p khf i+1 multiplied by some non- 
zero homogeneous 9 G k\Yo, . . . , Y n -i\ . Observe that the coefficients of rn p khj i+1 are 
homogeneous too. The computation of 9 and the remaining coefficients of X p k hf i+1 = 
9 ■ m p khj. +1 requires one more division free and well parallelizable straight-line program 

in k[Y ,..., Y n -i] of size . 

Putting everything together we obtain a division free straight-line program in 

k[Y ,...,Y n ^] 

of size Ldf^ + Aj and depth O(log20~i) + I + Aj which computes the homogeneous con- 
stant coefficient, say A G k[Yo, . . . , Y n —i] , of 9 ■ m p khj i+1 and the polynomial 9. We 

remark that A has degree at most (d5i)°^ . By Lemma || this coefficient has the form 
A = 9pf l ao for some m < ddi and is therefore nonzero. Since rrihj. +1 annihilates Vi+i 
modulo Jj we conclude that oo G ( h /j+i, Jj) holds in &[Yb> • • • , Y n ] = k[Xo, . . . , X n ] . 
Since the polynomial A can be evaluated by a division free straight-line program in 
k[Yo, . . . ,Y n -i\ of size L(dJj)°W + Aj and depth 0(log2d5i) + t + Aj , there exists by 
Lemma |6| a correct test sequence of length (Lb®^ + Aj) 3 in k n ~ l+1 for this complexity 
class. Because the polynomial A is nonzero we find in sequential time (L(d5i)°^ + Aj) 4 
and parallel time 0(log2d5i) + I + Aj a point 7 = (70, . . . , 7 n -i) G k n ~ l+l in this correct 
test sequence such that ^.(7) 7^ holds. Since .4 is homogeneous we may construct 
in the obvious way from the coordinates 70, • • • , 7n-i of 7 a linear change of the vari- 
ables Yq, . . . , Y n _i into new ones, say Yq, . . . , Y! n _ i , such that the homogeneous polyno- 
mial A' G fc[y ',...,l^_j] = /c[y ,---,^n-i] given by the equation A'(Yq, . . . = 
.4(Yo, • • • ) ^n-i) becomes monic in Y^_j . This implies that the homogeneous polynomial 
A (Yq, := a (Y , Y" n _j) is monic too in Y^_ { . 

Therefore, if we replace the variables Yq, . . . ,Y n by the new variables 

v' v' v' ■ v ■ -. v' ■ v 

2 0' • • • ' I n-ii 1 n-i+l - ~ J n-i+l> • • • j 1 n ■— 1 n 

we obtain a Noether normalization with respect to the ideal li+i = ( h fi+i, Jj) and hence 
with respect to Jj + i . 

In the next part of our algorithm we shall make use of the eliminating form ao instead 
of the "accidental" polynomial A which comes from our specific algorithm. We consider 
ao , A, 9 and pi as forms in the variables Yq, . . . , Y! n _ i which are related by the identity 



A_ 

Let p := 9p\ and observe that p is monic in Y^_j . The exponent k < d is known 
and the forms A , 9 and pi (and hence p = 9p\ ) are represented in our algorithm by a 
division free straight-line program £ in k\Y$, . . . , Y^_J = fc[Yo, . . . , Y n —i\ which has size 
L(d<5j)°W + Aj and depth 0{log2d5i)) + £ + Xi . The degrees of p and A are of order 
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Using N := max{deg[i, degA} = (do;) ^ parallel organized calls to the procedure E 
we interpolate A and ji with respect to the variable Yl a _ i in N many distinct points of 
k and compute their coefficients with respect to the variable (these coefficients 

are polynomials which belong to k[Y$, . . . , ). This can be achieved by a division 

free straight-line program in k\Y^ . . . , of size and depth (L + Ai)(d5i)°^ and 

0{log2d5i) + I + Xi respectively. Since the square matrix with entries in k\Y^ . . . , 
which correspond to the division (with remainder) of A by [i in the principal ideal domain 
k(YQ, . . . , Y^_j_ 1 )[Y^_j] is unimodular and since the identity oo = ^ holds, we are able to 
compute the coefficients of the polynomial ao with respect to Y! ri _ i (and hence ao itself) 
from the circuit E in sequential and parallel time (d5i)°^ and 0(log2d5i) respectively. 
Therefore we may suppose that ao (and its coefficients with respect to Y^ l _ i ) is given by 
a division free straight-line program in 1c[Yq, . . . , Y^_J = k[Yo, . . . , Y n _i] of size and depth 
(L + Aj)(d<5i)°W and 0{log 2 d5i) + I + \ respectively. 



5.1.2 Recursive Generation of a Primitive Element. 

The next recursive step is the computation of a primitive element Uj + i modulo Jj+i . We 
assume without loss of generality that the variables Yq, . . . ,Y n are already in Noether po- 
sition with respect to the ideal 2j+i = ( fi+i, Ji) (see "recursive Noether normalization"). 
So, we have the following integral ring extension 

R := k[Y , . . . , Y n ^} k[Y , Y n ]/( h f, l+1 , J, t ) = k[Y , Y n ]/(l i+1 ) 

Let us denote the fraction field of R by K (i.e. K := k(Yo, . . . , Y n -i~ %) ). It is clear 
that Xj + i generates a zero dimensional ideal in K[Y n _i, . . . , Y n ] . We are going to apply 



the method described in Section O. For this purpose we introduce new variables, say 
T n _j, . . . ,T n , and with respect to them for n — i < j < n the rings and fields Rj : = 
R[T n -U ■ ■ ■ , Tj-i, Tj+i, . . . , T n ] and Kj :— K(T n —i, . . . , Tj—i, ^j+i> • • • j T n ) introduced in 
[4.3| . For n — i < j < n we are going to construct polynomials hj E R[T] and gj G Rj[T] 
with the following properties : 

• hj £ R[T] is monic in T of degree at most d5{ in T and the polynomial hj(Yj) 
belongs to the ideal Ti+± 

• gj G Rj[T] is monic in T and has total degree at most d5i . The polynomial gj(Zj) 
belongs to the ideal Kj 0^ Zj + i 

Our aim is to represent the coefficients of hj and gj with respect to T by an evaluation 
scheme. 

With the notations and results of "recursive Noether normalization" there is already 
given an equation ao G = k[Yo, . . . ,Y n _i] for the integral dependence of Y n -i 

over R modulo the ideal 2j + i . This polynomial ao is given by a straight-line program in 
k® k R[Y n - i \ = k[Y ,...,Y n - i ] oflength (L + A i )(d5 i )° (1) and depth 0{log 2 (d5i) + £) + A; 
and it has degree d5i . Moreover it is homogeneous with respect to the natural grading of 

R[Y n -i] ■ 

Fix n — i < j < n and consider the variable Yj and the linear form Zj := T n ^iY n -i + 
\-Tj-iYj-i+Tj + iYj + i-\ \-T n Y n G Rj[Y n ^i, . . . ,Y n ] in the variables Y n _i, . . . , Y n (see 
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Section [4.3Q . To the homotheties given by the forms Yj and Zj in the Si -dimensional al- 
gebras K[Y n -i, . . . , Y n ]/(Ji) and Kj [Y n -i, . . . , Y n ]/(Ji) and to the bases of these algebras 
determined by l,Ui, . . . correspond matrices My } G K SixSi and Mz i G Kj lXSt . The 

characteristic polynomials of these matrices induce polynomials in the variables Y n -i,Yj 
and Y n _i,Zj respectively which have degree at most Si in these indeterminates. Since 
the variables Yq, . . . ,Y n are in Noether position with respect to Jj the coefficients of 
these polynomials belong to R and Rj respectively. These coefficients are homogeneous 
elements of degree at most Si of their respective graded rings. Let us denote the polyno- 
mials introduced in this way as Hj G R\Y n —i, Yj] and Gj G Rj\Y n —i, Zj] . They have the 
following properties : 

• Hj is monic in Yj and belongs to the ideal Jj 

• Gj is monic in Zj and belongs to the ideal Kj <S>_r Jj . 

These polynomials can be computed from the results of the circuit Tj as in "recursive 
Noether normalization" by means of the division free and well parallelizable algorithm || 
in sequential time S®^ . This means that their coefficients with respect to the variables 
Yj and Zj are given by division free straight -line programs in the rings k f^k R\Y n -i\ = 
k[Y , ... , Y n ^i] and k (g> k Rj [F„_j] = k[T n _i, ... , T j+ x, . . . ,T n , Y n _,] and that these 

straight-line programs have length 5f + Aj and depth 0(log2Si) + Aj . (Without going 
into details we remark here that one has to use the same type of subroutine to eliminate 
divisions by a certain polynomial which is monic in Y n _i as we did at the end of "recursive 
Noether normalization" when computing ao ) . 

For n — i<j<nwe are now going to construct the announced polynomials hj G R[T] 
and gj G Rj [T] . For this purpose we need the coefficients with respect to Y n -i of ao and 
the coefficients with respect to Y n _i , Yj and Y n ^i , Zj of Hj and Gj . This requires 
interpolation of these polynomials in dSi + 1 points of k which we have to substitute 
for the variable Y n _i (recall that Hj and Gj are already given by their coefficients with 
respect to Yj and Zj respectively). For this interpolation we need (2i + l)(dSi + 1) 
parallel organized calls to the whole procedure which originate a total sequential time cost 
of i(L + Ai)(dSi)°^ 1 ' and a parallel time cost of 0(log<idSi + €) + \ . 

Thus the coefficients with respect to the variable Y n _i of ao and for n — i < j < n 
the coefficients with respect to Y n _i , Yj and Y n _i , Zj of the polynomials Hj and Gj 
are represented by division free straight-line programs in k R = k[Yo, . . . , Y^,_j_i] and 
k <8>fc Rj = k[T n ^i, . . . , Tj-i, Tj+i, . . . , T n ] which have length i(L + Ai)(dSi)°^ and depth 
0(log2(dSi) + i) + Xi . Once these coefficients are given we are able to perform effectively 
computations in the graded k -algebras which for n — i < j < n are defined as follows : 

. Bj :=R[Y n ^,Yj]/(a ,Hj) 

• Cj := Rj{Y n _i, Zj]/ (ao, Gj) . 

Observe that Bj is a free R -module with basis represented by the set of monomials 
{Y^Yf :0<a<deg (a ), < < deg (Hj)} 
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and that Cj is a free Rj -module with basis represented by the set of monomials 
{Y^Zf :0<a<deg (a ), < < deg (Gj)}. 

The ranks of the free modules Bj and Cj are therefore bounded by d5f . Since the 
coefficients of ao with respect to Y n -i and the coefficients of Hj with respect to Y n -i , Yj 
are given, we know effectively (with no extra cost) the multiplication tensor of Bj . In the 
same sense the multiplication tensor of Cj is available. From this data we compute the 
characteristic polynomial Xyj £ R[T] of the homothety in Bj given by the multiplication 
by Yj and the characteristic polynomial Xzj G Rj[T] °f t ne homothety in Cj given by 
Zj . We put hj := Xy- an d 9j '■= Xz ■ One verifies immediately that hj and Qj have 
the required properties (for this purpose observe that ao £ , Hj £ Jj C Zj+i and 

G iiTj ® R Jj C g) B X i+ i holds). 

Applying Lemma and the observations at the end of Section 4T we obtain a linear 



form Ui + i which is a primitive element for the ideal generated by Xi+x in K\Y n -i, . . . , Y n ] , 
a homogeneous polynomial q £ R[T] and homogeneous polynomials p £ R and . . . , 

v n £ R[T] such that the following holds : q is monic and separable in T (hence square- 
free), q(v,i+i) belongs to the ideal yTj+i , p is nonzero and for n — i < j < n we have 
pYj — Vj(ui + {) £ y/Ti+i ■ Observe that we may suppose without loss of generality that p 
is monic in one of its variables, say Y n _j_i . 

These items (and their homogeneous coefficients in R = k[Yo, . . . ,Y n „j_i]) are rep- 
resented by a division free straight- line program in k[Yo, . . . ,Y n _j_i] whose length and 
depth are i(L + Ai){d8i)°^ and 0(log-2(iddi) + £) + Aj respectively. 



5.1.3 Clearing Extraneous Irreducible Components. 

In "recursive Noether normalization" we obtained a fairly explicit description of the projec- 
tive variety defined by the ideal li+i and of the localization of this ideal in K[Y n _i . . . , 1^,,] . 
However this projective variety may contain irreducible components at infinity (extrane- 
ous components) or components which count in Xj + i (or its localization) with higher 
multiplicities than one. 

The situation in the toric case is the same with "extraneous components at infinity" 

n 

replaced by "extraneous components contained in the union of hyperplanes V(J\ -Xi)"- 

%=\ 

We show now how extraneous components and multiplicities which appear in a natural 
way when cutting by the hypersurfaces V{ h f\), . . . , V( h f n ) , can be cleared out during the 
process. Let us concentrate upon the affine case. 

First recall that Ji+i is the intersection of all those homogeneous prime ideals of 
codimension i+1 which contain h fi+i and J« but do not contain X$. We observe also that 
the (homogeneous) primes of codimension i + 1 of R\Y n ^i, . . . , Y n ] = k\Yo, . . . , Y n ] which 
contain the ideal ( fx, ■ ■ ■ , h fi+i) but not Xq correspond to the primes in the localized 
ring K\Y n -i, . . . , Y n ] which contain the ideal generated by T%+\ but not Xq . We are going 
now to determine the irreducible components at finite distance of the projective variety 
defined by h fx, ■ ■ ■ , h fi+x (in other words defined by li+x ) doing computations in if . 
These components characterize geometrically the ideal Ji+x ■ For 1 < j < i + 1 let us 
write Fj := h fj(Y , . . . , Y„_i_i, . . . , ^f±^) £ K[u i+1 ] . Let h £ R[u i+1 ] be a 
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greatest common divisor of the polynomials Fi{ui+\), . . . , Fi + \{ui + \) and q(ui + ±) with 
respect to the principal domain . 

Since q(ui + i) is monic and separable in we may suppose without loss of generality 
that the polynomial h is monic and separable too in v,i+i . Moreover we may suppose 
that h is homogeneous with respect to the grading of . The coefficients of the 

polynomial h with respect to the variable Uj+i can be computed from the coefficients of 
Fi, . . . , Fi + i , q £ K[ui+i] by a well parallelizable and division free straight-line program 
of size (d5i)°W. This means that these coefficients can be computed by a division free 
straight-line program in k ®k R = k\Yo, . . . ,Y n -i] of size i(L + Ai)(d6i)°^ and depth 
0(log2(dSi)+£)+\i (here using the assumption that p is monic in Y^_j_i , we apply again 
the same trick as at the end of "recursive Noether normalization" in order to eliminate 
division by p ) . 

Observe now that the maximal ideals in K\Y n -i, . . . , Y n ] which contain the polynomial 
h(v,i+i) and for n — i < j < n the polynomials pYj — Vj(v,i+i) , are in one to one 
correspondence with the codimension i + 1 prime ideals in R[Y n -i, . . . , Y n ] which contain 
( h fi,..., h fi + i). Moreover h(ui + \) is squarefree. 

Note that Xo € k[Xo, . . . , X n ] = k[Yo, . . . , Y n \ is a linear form in Yo,...,Y n . Substi- 
tuting in this linear form for the variables y„_j, . . . , Y n the polynomials 

v n -i{u i+ i) v n (u i+1 ) 

, • • • , e a [ui+i\ 

p p 

and clearing the denominator p we obtain a representation in of the residue class 

of pXo modulo the ideal generated by Xi+i in K[Y n _i, . . . ,Y n ] . Let G G i?[itj + i] be this 
representation. Now we compute the greatest common divisor of G{ui + \) and h(ui + i) 
in the principal ideal domain A"[«j + i] . This greatest common divisor is represented by 
a squarefree polynomial h\ G which divides h G . Observe that hi is 

monic in Uj+i and homogeneous with respect to the grading of . The maximal 

ideals of K[Y n -i, . . . , Y n ] which contain and for n — i < j < n the polynomials 

pYj — Vj{ui + i) , are in one to one correspondence with the codimension i + 1 prime ideals 
of R[Y n -i, ...,Y n ] which contain (X , h f i+ i, . . . , h f n ) . 

Now it is easy to see that the maximal ideals in K[Y n -i, . . . ,Y n ] which contain for 
n — i < j < n the polynomials pYj — Vj{ui + \) and the homogeneous and monic in 
polynomial 

h(u i+1 ) 

q i+1 {u i+1 ) := — 

«i(%i) 

correspond exactly to the codimension i + 1 prime ideals which contain the radical ideal 
Jj+i (these are also the associated primes of Jj+i ). Observe that q^\ := q^xiT) is monic 
and separable with respect to the variable T . Finally the polynomials , . . . , Vn +1 ^ are 
defined (and computed) as the remainders of the division in R[T] of v n -i(T) . . . , v n (T) by 
the monic in T polynomial g , j + i(T) (note that this leaves the ideal in question unchanged). 
In this way we obtain an explicit description of an isomorphism K\Y n -i, . . . , Y n ]/ ( J«+i) — 
K[ui+i]/(qi+i(v,i+i)) = K[T]/(qi + i(T)) as wanted. Since hi is monic in Ui+i we are 
able to compute by the same trick as before qi + \ by a division free straight-line program. 
Because qi + \ is again monic in T the same argument applies to v^P , • • • , Vn +1 ^ . Fi- 
nally we observe that p and the coefficients with respect to the variable T of q{T) and 
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j (T) , . . . , vii (T) can be computed by a division free straight-line program Tin 
k\Yo, ■ ■ ■ , of size i(L + Ai)(d5i)°^ and depth 0(log2(d5i) + £) + Aj . Since c^+i 

is separable and monic with respect to T its discriminant (which we denote by pi+x ) is 
nonzero and belongs to R = k[Yo, . . . , . Since it is possible to compute pi + \ by a 

well parallelizable straight-line program in R of size (d5i)°^ we may suppose without 
loss of generality that the circuit T represents also the polynomial pi + % . 

We conclude this subsection with the remark that the toric case can be treated in 
exactly the same way replacing at the end of the construction for the affine case the linear 

n 

form Xq by the homogeneous polynomial T[ -^i ■ 

i=o 

5.2 Proof of Theorems 1, 2 and 3 

In this subsection we deduce from Lemma [l] and § and from Proposition [l| first Theorem 
H and then Theorem [l] and ||. 

Proof of Theorem pj. Let assumptions and notations be the same as in the statement 
of Theorem ^ and as at the beginning of this section. The algorithm underlying the proof 
proceeds in n recursive steps. 

Let us now show Theorem [| {%) which deals with the affine case. In the first step of 
our algorithm we apply Proposition [l] just to the input straight -line program (5 and the 
input polynomials fi, ■ ■ ■ , f n in the following sense : we put formally i = and assume 
To to be the empty circuit. 

Suppose now 1 < i < n and let us consider the [i + l)-th step of our procedure. Let 
be given an arithmetic network M% with parameters in k which produces the input for the 
(i + l)-th recursion step of our procedure, namely a division free straight-line program 
Pi in k[YQ, . . . ,Y n ^i] which evaluates pi and the coefficients with respect to T of the 
polynomials qi , v^_ i+1 , . . . , Vn . Suppose also that 0i contains the information about 
the coefficients of the linear form Ui which for example may be stored as parameters of 
Pi . Assume inductively that the polynomial qi is monic and separable with respect to 
the variable T and that pi is its discriminant. Let Li and li be the nonscalar size and 
depth of Mi respectively. Note that the size and depth of the circuit Pi do not exceed Li 
and li respectively. 

To the circuits P and Pi, the polynomials Pi , qi, v^-i+l! • • • > w ^ an< ^ 

the A; -linear form itj we apply the compression algorithm underlying Lemma [lj. We 
obtain a new division free straight-line program Ti in k[Yo, . . . ,Y n —i] which represents 
the coefficients of u% , the polynomial pi and the coefficients with respect to T of qi 

(i) (i) n 

and v n _ i+l , . . . ,Vn ■ From Lemma |lj we deduce that the circuit Ti has size and depth 
0((i 5 + L)5} 1 ) and 0((log 2 i + t)log 2 8i) . 

We apply now to the circuit P and Tj , the polynomials 

fl, • • • , fi+i, Ph Qi, «i-i+u • • • , 4 l) 

and the k -linear form U{ the elimination procedure underlying Proposition [l|. We obtain 
a linear change of the variables Yq, . . . ,Y n of the i— th step into new variables Yq, . . . ,Y^ 
such that Yq, . . . , Y^ are in Noether position with respect to the ideal Jj+i , the vari- 
ables Yq, . . . ,Y^_i_i being free. Moreover we obtain a new k -linear form Uj+i and new 
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polynomials p i+1 G k[Y^, . . . , F^^J and % + i , , • • • , v { n +l) G fc[F ', . . . , F^_j_ l5 T] . 

All these items are represented by a division free straight-line program in k\Y$, . . . , 
• • • , which is generated by an arithmetic network Aj+i with parameters in k . 

From Proposition H we deduce that the size and depth of are 0(i(L+(i 5 +L)5j 1 )(d5i)°^) = 
i 6 L(d5 t )°^ and 0((log 2 (d5i) + £) + (log 2 i + ^Zo^) = 0(log 2 (id) + £)% 2 ^) . 

Similarly the size and the depth £i+\ of Aj are (iLd5i)°^ +L{ and 0(log 2 (id) + 
l)log 2 8i) + 4 . 

Finally we obtain for i := n — 1 the output of the algorithm underlying Theorem |3|(i) . 
This output consists in the coefficients of a nonzero linear form u G k[X\, . . . , X n ] and of 
one-variate polynomials q , v\, . . . ,v n G k[u] such that the following holds : 

• degq = 5 n 

• max{degvi, . . . , degv n } < 5 n 

• - ..,/n) = (<l, x l ~ Vi,...,,X n - v n ). 

The output is represented by the arithmetic network N n whose parameters belong to 
k and whose size and depth are (nLd5)°^ and 0(n(log 2 (nd) + £)log 2 5) respectively. 

The proof of Theorem || (ii) which deals with the toric case is completely analogous. 
One has just to replace in it the application of Lemma [2] by the application of Lemma |l| 
in the same sense. In order to avoid repetitive arguments we omit this proof. ■ 

^,From Theorem ^ (i) and (ii) we deduce now easily Theorem [l] and ||. 

Proof of Theorem [J and [|. Let us first concentrate upon the afflne case, namely 
Theorem |l[ 

Let assumptions and notations be the same as in the statement of Theorem [l]. Applying 
Theorem || (i) to this situation we obtain an arithmetic network M with parameters in k 
which computes the coefficients of a & -linear form u and of polynomials q , v% , . . . , v n G 
k[T] satisfying the following conditions : 

(*) degq = 5 n and max{degvi ; I < i < n} < 5 n 
(**) (/i,...,/„) = (q(u),Xi - Vl (u),...,X n -v n (u)) 

The network M has size (nd5L)°^ and depth 0(n(log 2 (nd) + £)log 2 5) . 

Let H = ct\Xi + • • • + a n X n be the representation of the nonzero linear form H G 
k[Xi, . . . , X n ] by its coefficients a%, . . . , a n G k and let p be the characteristic (or min- 
imal) polynomial of the A: -linear endomorphism of k[T]/(q) induced by the polynomial 

aivi H \-a n v n G k[T] . We can compute the coefficients of p from the data c^i, . . . , ct n , 

q and v\, . . . ,v n by a well parallelizable algorithm in sequential time 5® ■ Thus we may 
suppose without loss of generality that M computes also the coefficients of the monic poly- 
nomial p. The condition (**) above expresses an isomorphism between the A; -algebras 
k[T]/(q) and k[X\, . . . , X n ]/ . . . , f n ) . From the particular form of this isomorphism 
we infer that p is also the characteristic (or minimal) polynomial of the k -linear endo- 
morphism of k[Xi, . . . , X n ]/(fi, . . . , f n ) induced by the multiplication by H . This implies 
that p(H) G (/1, . . . , f n ) holds, whence the conclusion of Theorem [[]. 

Theorem || is deduced from Theorem ||| (ii) in almost textually the same manner. ■ 
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Lemma [T] and Proposition [l] can be combined in the same way as in the proof of 
Theorem || {%) in order to obtain the following elimination result for reduced complete 
intersection ideals of positive dimension : 

Proposition 2 Let fi,...,f s be polynomials which belong to k[X±, . . . , X n ] . Suppose 
that fi,...,f s form a regular sequence in k[X\, . . . , X n ] . For any 1 < j < s let 
Sj := degV(fi, . . . , fj) be the geometric degree of the affine variety defined by the ideal 
(/i, . . . , fj) which we assume to be radical. Write r := n — s , 5 := max{5j ; 1 < j < s} 
and d := max{degfj ; 1 < j < s} . Suppose that the polynomials fx,...,f s are given 
by a division free straight-line program (5 in k[X±, . . . , X n ] of length L and depth i. 
Then there exists an arithmetic network with parameters in k which has size (sd5L) ^ 
and depth 0(s(log2{sd) + t)log 2 S) and which produces from the circuit as input, the 
following items : 

• a nonsingular matrix of k nxn which is given by its coefficients and which transforms 
the variables X\, . . . , X n into new ones, say Y±, . . . , Y n , 

• a nonzero k -linear form u £ k[Y r+ \, . . . , Y n ] , 

• a division free straight-line program 7 in k[Yi, . . . ,Y r ] which represents a nonzero 
polynomial p £ k[Yi, . . . , Y r ] and the coefficients with respect to the variable u of 
certain polynomials q, v\, . . . , v n £ k[Y\, . . . , Y r , u] . 

These items have the following properties : 

i) the new variables Y\,...,Y n are in Noether position with respect to the ideal 
. . . , f s ) , the variables Yi,...,Y r being free, 

ii) the polynomial q is monic and separable in u and p is its discriminant, 
Hi) in the localized ring k[X\, . . . ,X n ] p the following ideals are identical 

(fi,---,fs) P = {q,pXi -vi,...,pX n -v n ) p , 

iv) the polynomials q, p and v\,...,v n satisfy the degree conditions degq = deg u q = 
5 S , degp < 5^ , max{deg u Vj ; 1 < j < n} < 5 S and max{degvj ; 1 < j < n} < 25$ , 

v) the length of the straight-line program 7 is 0{{s b + L)5\ l ) = (sL5)°W and its depth 
is 0((log 2 s + £)log2S s ) = 0((log 2 s + i)log 2 5) ■ 

We omit the proof of this proposition which is essentially the same as the proof of 
Theorem ||] (i) . 

Of course Proposition [2] admits also a toric version with almost the same proof as 
Theorem |3| (ii) . We are going now to formulate a slightly more general statement in which 

n 

the role of the form Yl Xi is played by an arbitrary nonzero polynomial g £ k [X\ , . . . , X n ] . 

We think this statement is interesting by its own because it says that our algorithms are 
able to "avoid" undesired points which are contained in a previously given hypersurface 
V{g). 
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Let be given nonzero polynomials fx, . . . , f s ,g S k[Xx, . . . , X n ] with s < n. We say 
that f\, . . . ,f s form a secant family ("suite secante") which avoids the hypersurface V(g) 
if for any 1 < j • < s the localization V(fx, ■ ■ ■ , fj) g is not empty and if for any irreducible 
component C of V(fx, ■ ■ ■ , fj) the following implication holds : 

if C is not contained in the hypersurface V(g) then dimC = n — j . 

Let 1 < j < s . We say that an irreducible component C of V(/i, • • • , fj) avoids the 
hypersurface V(g) if C is not contained in V(g) . The union of all irreducible compo- 
nents of V(fx, ■ ■ ■ , fj) which avoid the hypersurface V(g) is called the part of the affine 
variety V(/i, • • • , fj) which avoids the hypersurface V(g) and we denote this part by 
V'(fx, ■ ■ ■ , fj) ■ We write deg'V(fx, ■ ■ ■ , fj) for its geometric degree. Thus V'(fx, ■ ■ ■ , fj) 
is the Zariski closure of the locally closed subset V(/i, • • • , fj) g of A n and we have 
deg'V(fx, ■ ■ ■ , fj) = degV'(fx, ■ ■ ■ , fj) ■ Moreover if fx, ■ ■ ■ , f s form a secant family which 
avoids the hypersurface V(g) then deg'V(fx, ■ ■ ■ , fj) > holds for any 1 < j < s . 

With these notions and notations we are now able to formulate our statement which 
generalizes Theorem ^ (ii) to positive dimension : 

Proposition 3 Let fx, ■ ■ ■ , fs and g be polynomials of k[Xx, . . . , X n ] such that fx, ■ ■ ■ , fs 
form a secant family which avoids the hypersurface V(g) ■ For any 1 < j < s let 5'j := 
deg'V(fx, • • • , fj) be the degree of the part of the affine variety V(fx, ■ ■ ■ , fj) which avoids 
the hypersurface V(g) and suppose that the localized ideal (fx, ■ ■ ■ ,fj)g is radical. Write 
r := n — s , 5' := max{5j ; 1 < j < s} and d := max{degfj ; 1 < j < s} . Assume 
that the polynomials fx, ■ ■ ■ , f s and g are given by a division free straight-line program (3 
in k\Xx, ■ ■ ■ ,X n ] of length L and depth £ . Then there exists an arithmetic network with 
parameters in k which has size (sd5' L)°^ and depth 0(s(logsd + £)log25') and which 
produces from the circuit (3 as input, the following items : 

• a nonsingular matrix of k nxn which is given by its coefficients and which transforms 
the variables Xx, ■ ■ ■ , X n into new ones, say Yx , . . . , Y n , 

• a nonzero k -linear form u G k[Y r+ x, ■ ■ ■ , Y n ] , 

• a division free straight-line program 7 in k[Yx, . . . ,Y r ,u] which represents a nonzero 
polynomial p G k[Yx, . . . ,Y r ] and the coefficients with respect to u of certain poly- 
nomials q, vx , . . . , v n G k [Yx , . . . , Y r , u] . 

These items have the following properties : 

i) the new variables Yx, . . . ,Y n are in Noether position with respect to V'(fx, ■ ■ ■ , fs) > 
the variables Yx,...,Y r being free, 

ii) the polynomial q is monic and separable in u and p is its discriminant, 
Hi) in the localized ring k[Xx, ■ ■ ■ ,X n ] p . g the following ideals are identical 

(fl, ■ ■ ■,fs)pg = (q,pXx - Vx,-. -,pX n - V n ) p . g , 

iv) the polynomials q, p and vi,...,v n satisfy the degree conditions degq = deg u q = 
5' s , degp < 2<5 S 3 , max{deg u Vj ; 1 < j < n} < 5' s and max{degvj ; 1 < j < n} < 
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v) the length of the straight-line program 7 is 0((s 5 + L)5' s 11 ) = (sL5')°^ and its 
depth is 0((log 2 s + l)log 2 8' s ) = 0((log 2 s + l)log 2 8') . 

We omit the proof of this proposition which essentially the same as the proof of The- 

n 

orem |3] (ii) with the form FJ Xi replaced by the polynomial g . 

i=l 

6 Division Theorems 

In this section we explain how Proposition ^ can be applied in order to obtain refined 
complexity and degree bounds in the division theorems of [19|, |^(J an d |52|] (see also [26] 
and |41[| ). Let us first state our results. 

Theorem 5 Let f\, . . . ,f s and g be polynomials of k[X\, . . . , X n ] such that fi, ■ ■ ■ , f s 
form a regular sequence and g belongs to the ideal (/1, . . . , f s ) . For 1 < j < n denote 
by 5j := degV(f\, . . . , fj) the geometric degree of the affine variety defined by the ideal 
(/l,... , fj) which we suppose to be radical. Write 5 := max{5j ; 1 < j < s — 1} and 
d:= max{degfj ; 1 < j < n} . Assume that the polynomials fi, . . . , f s and g are given by 
a division free straight-line program (3 in k[X±, . . . ,X n ] of length L and depth i. Then 
there exists an arithmetic network with parameters in k which has size (sd5L)°W and 
depth 0(s(log 2 sd + £)log 2 5) and which produces from the circuit (3 as input a division 
free straight-line program 7 in k[X\, . . . ,X n ] such that 7 represents certain polynomials 
pi,...,p s G k[X%, . . . , X n ] with the following properties : 



i) 9 = Pl/l H ^Psfs 

ii) max{degpj ; 1 < j < s} < (2s 2 d + max{d, degg})5 

Moreover the circuit 7 has size (degg) 2 (sd5L)°^ and depth 0(s(log 2 sd + £)log 2 5) . 
Let us remark that allowing divisions in the circuit 7 we may diminuish its size to 
(sd5L)°W . 



Let fx, . . . , f s be nonconstant polynomials representing an equation system of k[X\, 
. . . , X n ] . Define polynomials g%, . . . ,g n £ k[Xi, . . . , X n ] as follows : if the characteristic of 
k is zero, choose g\, . . . ,g n as generic linear combinations of the polynomials fi, ■ ■ ■ , f s 
and if the characteristic is positive, choose gx,...,g n as generic linear combinations of 
the set of polynomials {Xifj ; 1 < i < n , l<i<s}- For 1 < i < n denote by 
5i := degV(g\, . . . , g n ) the geometric degree of the affine variety V(g±, . . . , gi) . Then we 
call S:= max{5i ; 1 < i < n} the generic geometric degree of the equation system given by 
/1, . . . , f s . Observe that the generic geometric degree of the system fi,...,f s is positive 
even if fi,...,f s generate the trivial ideal. 

With these notion we can state our next result, namely 

Theorem 6 (effective Nullstellensatz) Let fi, . . . , f s be polynomials of the ring 
k[X\, . . . , X n ] which generate the trivial ideal. Let 5 be the generic geometric degree of the 
equation system fi, . . . , f s and let d := max{degfj ; 1 < j < s} . Suppose that the poly- 
nomials /1, . . . ,/ s are given by a division free straight-line program (3 in k[X\, . . . ,X n ] 
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of length L and depth £ . Then there exists an arithmetic network with parameters in k 
which has size (ndSL) ^ 1 ' and depth 0(n(lognd + £)log2S) and which produces from the 
circuit (5 as input a division free straight-line program 7 in k [X\ , . . . , X n ] such that 7 
represent certain polynomials pi, . . . ,p s £ k[X±, . . . , X n ] which have the following proper- 
ties : 



i) 1 =Pi/i H \-Psfs 

ii) max{degpj ; 1 < j < s} < 3n 2 dd 

Moreover the circuit 7 has size (nd5L)°^ and depth 0(n(lognd + £)log2S) . 

If we replace in the statement of Theorem || and || the parameter 5 by the Bezout 
estimations d s and d n respectively we obtain the worst case complexity bounds of [19], 



Theoreme 5.1 and 5.2, J4(| Theorem 1 and 2 and |y], Lemma 15 and 16, Theorem 19 and 
24. 

Let us now make some comments first on the degree and then on the complexity 
bounds of Theorem || and |6] and their proofs. The degree bounds in Theorem [5] (ii) and 
^ (ii) are noteworth by their own. They are due to the paper |5^] from which they follow 



immediately by performing the following modification in the proof of [52], Lemma 15 and 
16 : 



every time when the Bezout inequality ([29], Theorem 1) is applied in order to produce 
a degree bound of d? , this bound is replaced by the number 5j , where 1 < j < s or 
1 < j < n, following the context (see [53]). 

We are now going to explain how the complexity bounds of Theorem |5| (i) and |6] (i) 
follow from Proposition |2|. In order to make explanations easier we refer only to proofs in 
paper [19] and to the sequential complexity model (compare with p0[ for considerations 
of nonscalar parallel complexity). We obtain a proof of Theorem || (i) and || (i) almost 
directly from the proof of |ll| Theoreme 5.1. and 5.2. if we perform systematically the 
following modifications in the paper [|l9| : Proposition 2.3.1. and Proposition principale 
2.4.1. are replaced everywhere they are used in [19] by our Proposition]^. Moreover in all 



degree estimations, when the application of the Bezout inequality produces a d? bound 
(for 1 < j < s or 1 < j < n), this bound is replaced by 5j . With these modifications the 
proof of our Theorem || (i) and ^ (i) is almost textually the same as that of | |19| ] , Theoreme 
5.1 and 5.2. There is only one point where caution is necessary : the application of duality 
techniques in [|H]] is based on the simple minded decomposition formula for duality (3) 
in Section 3.2 of the same paper. The size of this decomposition may happen to be too 
big for our purpose here. In order to remedy this problem we have to use a different view 



of duality which is close to [52 and |53|| . This will be contained in a forthcoming paper 
where complete and self contained proofs of Theorem || and || will be given. 

For getting simultaneously with the complexity bounds also the degree estimations in 
Theorem |5| and |6] one has to adapt the algorithm underlying the proof of Theorem || (i) 
and |6| (i) to the particular constructions in 52] and ]53[ ]. 
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7 Conclusions 



This final section is devoted to the announcement of further complexity results which can 
be obtained by our method. We hope that these results, which are intermediary products 
of work still in progress, shed light on the possibilities and limitations of our approach. 
First let us turn back to our discussion at the end of Section |3[ 

Remark || and Proposition [l] imply that our main results, namely Theorem |], ||, ||L ||, ^ 
and Proposition ||, || remain true if we replace in their conclusions "arithmetic network with 
parameters in k " by "arithmetic network over k with factorization gates" and "straight- 
line program in k[Xi, . . . ,X n ]" by "straight-line program in k[Xi, . . . , X n ] ". Parallel 
to this observation we may ask what happens to our algorithms when we try to transfer 
them to the bit complexity model (here we suppose that A; is Q or a finitely generated 
transcendental extension of Q). We announce here that our procedures allow perfect a 
modularization modulo suitably (randomly) chosen primes of not too big height. This 
implies that a suitably randomized version of our main complexity results remain valid 
in the bit model if we take into account the (bit) size of the rational parameters which 
appear in the input circuit . Finally let us state just one application of our method to 
an elimination problem in semialgebraic geometry. 

Let k be an ordered hilbertian field with factorization at moderate costs and let K 
be a real closure of k. Let k:= K := K(i) with i 2 = —1 be the corresponding algebraic 
closure of k and K . We consider the affine space A n := A n (k) equipped with the Zariski 
topology whose closed sets are the &; -definable algebraic subsets of k n . Let W C A n be 
a closed subset and let W = C\ U • • • U C s be its decomposition in irreducible components 
with respect to this topology. Thus W and C\, . . . ,C S are k— definable algebraic subsets 
of k n . Let 1 < j < s and consider the irreducible component Cj of W . We call Cj a 
real component of W if the real variety Cj D K n contains a smooth point of Cj . Let 

I = {j E IN : 1 < j < s, where Cj is a real component of W}. 

We call the affine variety W" := U Cj the real part of W and we define the real degree 

jei 

of Was deg"W:= degW" = £ degCj . Observe that deg"W = holds if and only if the 

jei 

real part W" of W is empty. 

Let / be a nonconstant and squarefree polynomial of k[X±, . . . , X n ] which defines 
a hypersurface W := V{f) of A n . Let V := W D K n the real variety given by the 
polynomial / . Suppose that V is nonempty, bounded and smooth with regular equation 
f . Assume that the variables X\, . . . ,X n are in generic position with respect to V and 
consider for any < i < n the polar variety of W corresponding to the linear space 
defined by the equations X±, . . . } Xi . Denote by Wi this polar variety and consider the 
real variety Vi := Wi n K n . With these notions and notations it is not too difficult to 
deduce from the Weak Transversality Theorem (a la Sard-Thorn) that for any < i < n 
the following facts are true : 

• W = W , 

• W{ is a nonempty equidimensional affine subvariety of A n which is smooth in all 
of its points which are smooth with respect to W , 
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• the real part W" of the polar variety Wj coincides with the Zariski closure of Vi in 



Let us write 6" := deg"Wi for < i < n and 8" := max{8'{ ; < % < n} and 
d:= degf . 

With these notations and assumptions we have the following complexity result : 

Theorem 7 (in collaboration with B. Bank and R. Mandel) Suppose that the equation 
f is given by a division free straight-line program (3 in k [X\ , . . . , X n ] of length L and 
depth £ . Then there exists an arithmetic network in k with factorization gates which 
has size (nd8" L)°^ and depth 0(n(log2(nd) + l)log28") which from the circuit [3 as 
input produces the coefficients of a nonzero linear form u G k[X\, . . . , X n ] and of certain 
polynomials q, v i, . . . , v n 6 k[u] . These polynomials have the following properties : 

i) q is monic and squarefree and has degree 8" . Moreover v±, . . . ,v n satisfy the degree 
bound max{degvj ; 1 < j < n} < 8'^ 

ii) For each semialgebraically connected component C of V there exists a point £ G C 
and a field element r G K such that q(r) = and £ = («i(r), . . . ,v n (r)) holds. 

We observe here that the network computes only points £ which are critical points of the 
projection map of V into K induced by the coordinate X n . 

Acknowledgments : We thank the unknown referee for his inspiring comments and 
his suggestions to improve the presentation of this paper. 

References 

[1] M.E. Alonso, E. Becker, M.-F. Roy, T. Wormann : Zeroes, multiplicities and idempotents for zerodi- 
mensional systems, to appear in Proc. MEGA '94, Birkhauser Progress in Mathematics. 

[2] W. Baur, V. Strassen : The complexity of partial derivatives, Theoret. Comput. Sci. 11 (1982) 
317-330. 

[3] C. Berenstein and A. Yger : Effective Bezout identities in Q[Xi, . . . , X n ] , Acta Math., 166 (1991) 
69-120. 

[4] C. Berenstein and A. Yger : Une formule de Jacobi et ses consequences, Ann Sci. E.N.S., J > %eme 
serie, 24 (1991) 363-377. 



• Vi is defined by the equations f, ■ ■ ■ , q^- , 

• for any i < j < n the ideal 




is radical. 



42 



[5] S.J. Berkowitz : On computing the determinant in small parallel time using a small number of 
processors, Inf. Proc. Letters 18 (1984) 147-150. 

[6] J. Briangon : Sur le degre des relations entre polynomes, C. R. Acad. Set. Paris, Serie I Math. 297 
(1982) 553-556. 

[7] D. W. Brownawell : Bounds for the degree in the Nullstellensatz, Annals of Math. 126 (1987) 
577-591. 

[8] L.Caniglia, A.Galligo and J.Heintz : Some new effectivity bounds in computational geometry, Proc. 
AAECC-6, T. Mora, ed., Springer LNCS 357(1989) 131-152. 

[9] L. Caniglia, A. Galligo and J. Heintz : Borne simplement exponentielle pour les degres dans le 
theoreme des zeros sur un corps de caracteristique quelconque, C.R. Acad. Sci. Paris, Serie I, Math. 
307 (1988) 255-258. 

[10] J. Canny : Some algebraic and geometric problems in PSPACE, Proc. 20-th Ann. ACM Symp. Theory 
of Computing (1988) 460-467. 

[11] A.L. Chistov : Polynomial -time computation of the dimension of components of algebraic varieties 
in zero-characteristic, Preprint Universite Paris Val de Marne. 

[12] A.L. Chistov, D.Yu. Grigor'ev : Subexponential time solving systems of algebraic equations, LOMI 
Preprints E-9-83, E- 10-83, Leningrad (1983). 

[13] F. Cucker, M. Karpinski, P. Koiran, T. Lickteig, K. Werther : On real Turing machines that toss 
coins, manuscript University Bonn (1995) 

[14] P. Debes : On the irreducibility of the polynomials P(t m ,Y), Journal of Number Theory, vol. 42 (2) 
(1992) 141-157 

[15] R.A. De Millo, R.J. Lipton : A probabilistic remark on algebraic program testing, Inf. Proc. Letters, 
vol. 7 (4) (1978) 193-195. 

[16] A. Dickenstein, N. Fitchas, M. Giusti, C. Sessa : The membership problem of unmixed ideals is 
solvable in single exponential time, Discrete Applied Mathematics 33 (1991) 73-94. 

[17] T.W. Dube : A combinatorial proof of effective Nullstellensatz, J. Symb. Comp. 15 (1993) 277-296. 

[18] D. Duval : Evaluation Dynamique et cloture algebrique en Axiom, Journal of Pure and Applied 
Algebra 99 (1995) 267-295. 

[19] N. Fitchas, M. Giusti and F. Smietanski : Sur la complexity du theoreme des zeros, in Approxima- 
tion and Optimization in the Caribbean II, Proceedings of the Second International Conference on 
Approximation and Optimization, La Habana, 1993, M. Florenzano et al., eds., Approximation and 
Optimization, Peter Lang Verlag, Frankfurt (1995) 274-329. 

[20] M. Fried : On the Sprindzuk-Weissauer approach to universal Hilbert subsets, Israel J. Math. 51 
(4) (1985) 347-363. 

[21] W. Fulton : Intersection Theory, Ergebnisse der Mathematik, 3 Folge Band 2, Springer Verlag (1984). 



43 



[22] C.B. Garcia and W. I. Zangwill : Pathways to solutions, fixed points and equilibria, Prentice-Hall, 
N.J. (1981) 

[23] J. von zur Gathen : Parallel arithmetic computations : a survey, Proc. 13-th Conf. MFCS. Springer 
LNCS 233(1986) 93-112. 

[24] M.Giusti and J.Heintz : La determination des points isoles et de la dimension d'une variete algebrique 
peut se faire en temps polynomial, In Computational Algebraic Geometry and Commutative Algebra, 
Proceedings of the Cortona Conference on Computational Algebraic Geometry and Commutative 
Algebra, D. Eisenbud and L. Robbiano, eds., Symposia Matematica, vol. XXXIV, Istituto Nazionale 
di Alta Matematica, Cambridge University Press (1993). 

[25] M. Giusti, J. Heintz, J.E. Morais, L.M. Pardo : When polynomial equation systems can be "solved" 
fast? in Proc. 11th International Symposium Applied Algebra, Algebraic Algorithms and Error- 
Correcting Codes, AAECC-11, Paris 1995, G. Cohen, M.Giusti and T. Mora, eds., Springer LNCS 
948 (1995) 205-231. 

[26] M. Giusti, J. Heintz and J. Sabia : On the efficiency of effective Nullstellensatze, Computational 
Complexity 3 (1993) 56-95. 

[27] L. Gonzalez Vega, G. Trujillo : Using symmetric functions to describe the solution of a zero- 
dimensional ideal, in Proc. 11th International Symposium Applied Algebra, Algebraic Algorithms and 
Error-Correcting Codes, AAECC-11, Paris 1995, G. Cohen, M. Giusti and T. Mora, eds., Springer 
LNCS 948 (1995) 232-247. 

[28] K. Haegele, J.L. Montana : Polynomial random test for the equivalence problem of integers given by 
arithmetic circuits. Manuscript Universidad de Cantabria (1995). 

[29] J. Heintz : Fast quantifier elimination over algebraically closed fields, Theoret. Comp. Sci. 24 (1983) 
239-277. 

[30] J. Heintz : On the computational complexity of polynomials and bilinear mappings, in Proc. 5th Inter- 
national Conference Applied Algebra, Algebraic Algorithms and Error-Correcting Codes AAECC-5, 
Menorca 1987, L. Huguet and A. Poli, eds., Springer LNCS 356 (1989) 269-300. 

[31] J. Heintz and J. Morgenstern : On the intrinsic complexity of elimination theory, J. of Complexity 9 
(1993) 471-498. 

[32] J. Heintz and CP. Schnorr : Testing Polynomials which are easy to compute, Proc. 12th Ann. ACM 
Symp. on Computing (1980) 262-268; also in Logic and Algorithmic. An International Symposium 
held in Honour of Ernst Speaker, Monographie No. 30 de l'Enseignement de Mathematiques, Geneve 
(1982) 237-254. 

[33] J. Heintz and M. Sieveking : Absolute primality of polynomials is decidable in random polynomial 
time in the number of the variables, 8th International Colloquium on Automata, Languages and 
Programming ICALP 81, Springer LNCS 115 (1981) 16-28. 

[34] O. Ibarra, H. Moran : Probabilistic algorithms for deciding equivalence of straight-line programs, J. 
ACM 30, Vol. 1, 217-228 



44 



[35] B. Iversen. Generic Local Structures of the Morphisms in Commutative Algebra. Springer Lect. Notes 
in Math. 310, Springer- Verlag (1973). 

[36] J. Ja' Ja' : An introduction to parallel algorithms, Addison-Wesley (1992) 

[37] E. Kaltofen : Greatest common divisors of polynomials given by straight line programs, J. ACM 35 
(1) (1988) 234-264. 

[38] E. Kaltofen : Factorization of polynomials given by straight-line programs, Randomness in Compu- 
tation, Advances in Computing Research 5, S. Micali, ed., JAI Press Inc., Greenwich, CT. (1989) 
375-412. 

[39] J. Kollar : Sharp Effective Nullstllensatz, J. of A. M.S. 1 (1988) 963-975. 

[40] T.Krick and L.M.Pardo : A Computational Method for Diophantine Approximation, to appear in 
Proc. MEGA '94, Birkhauser Progress in Mathematics. 

[41] T. Krick, L.M. Pardo : Une approche informatique pour l'approximation diophantienne, C. R. Acad. 
Set. Paris, t. 318, Serie I, no. 5 , (1994) 407-412. 

[42] E. Kunz : Kdhler Differentials, Advanced Lectures in Mathematics, Vieweg Verlag, Braun- 
schweig/Wiesbaden (1986) 

[43] Y.N. Lakshman, D. Lazard : On the complexity of zero-dimensional algebraic systems, Proc. 
MEGA '90, Carlo Traverso, ed., Prog. Math., vol. 94, Birkhauser Boston (1991) 217-225. 

[44] D. Lazard : Resolution des systemes d'equations algebriques, Theor. Comp. Sci. 15 (1981) 77-110. 

[45] A.K. Lenstra, H.W. Lenstra, L. Lovasz : Factoring polynomials with rational coefficients, Math. Ann. 
261 (1982) 515-534. 

[46] F.S. Macaulay : The Algebraic Theory of Modular Systems, Cambridge University Press (1916) 

[47] M. Moeller : Systems of algebraic equations solved by means of endomorphisms, in Proc. 10th Inter- 
national Conference Applied Algebra, Algebraic Algorithms and Error-Correcting Codes AAECC-10, 
Puerto Rico 1993, G. Cohen, T. Mora and O. Moreno, eds., Springer LNCS. 

[48] J.L. Montana and L.M. Pardo : Lower Bounds for Arithmetic Networks, AAECC, vol. 4 (1993), 
1-24. 

[49] J. Morgenstern : How to compute fast a function and all its derivatives, Prepublication No. 49, 
Universite de Nice 1984. 

[50] L.M. Pardo : How lower and upper complexity bounds meet in elimination theory, to appear in Proc. 
11th International Symposium Applied Algebra, Algebraic Algorithms and Error-Correcting Codes, 
AAECC-11, Parisl995, Springer LNCS 948 (1995) 33-69. 

[51] P. Philippon : Denominateurs dans le theoreme des zeros de Hilbert, Acta. Arith. 58 (1991) 1-25. 

[52] J. Sabia and P. Solerno : Bounds for traces in complete intersections and degrees in the Nullstellen- 
satz, in AAECC, vol. 6 (1995) 353-376. 



45 



[53] J. Sabia, P. Solerno : On intrinsic bounds in the Nullstellensatz, Manuscript Univ. de Buenos Aires 
(1995). 

[54] A. Schonhage : On the power of random acces machines, in Proc. 6th International Conference on 
Automata, Languages and Programming, ICALP, Graz (1979) 520-529 

[55] J.T. Schwartz : Fast probabilistic algorithms for verification of polynomial identities, J. ACM 27 
(1980) 701-717 

[56] M. Shub, S. Smale : Complexity of Bezout's theorem V : Polynomial time, Theoretical Comp. Sci. 
133 (1994) 141-164. 

[57] M. Shub, S. Smale : On the intractability of Hilbert's Nullstellensatz and an algebraic version of NP 
/ PI, IBM Research Report, Yorktown Heights (1994) 

[58] V.G. Sprindzuk : Diophantine equations involving unknown primes, Trudy MIANSSR 158 (1981) 
180-196. 

[59] V. Strassen : Algebraic Complexity Theory, Handbook of Theoretical Computer Science, ch. 11 (1990) 
634-671. 

[60] V. Strassen : Vermeidung von Divisionen, Crelle J. Reme Angew. Math. 264 (1973) 184-202. 

[61] V. Strassen : Die Berechungskomplexitat von elementarsymmetrischen Funktionen und von Interpo- 
lationskofnzienten, Numer. Math. 20 (1973) 238-251. 

[62] M. Yasumoto : Nonstandard arithmetic and Hilbert's irreducible sequences, manuscript (199?). 

[63] P. Zariski, O. Samuel : Commutative Algebra II. 

[64] R.E. Zippel : Probabilistic Algorithms for Sparse Polynomials, Springer LNCS 72 (1979) 216-226. 



46 



